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ABSTRACT 

We present a new multi-component spectral energy distribution (SED) decomposition method and use it 
to analyze the ultraviolet to millimeter wavelength SEDs of a sample of dusty infrared-luminous galaxies. 
SEDs are constructed from spectroscopic and photometric data obtained with the Spitzer Space Telescope in 
conjunction with photometry from the literature. Each SED is decomposed into emission from populations 
of stars, an AGN accretion disk, PAHs, atomic and molecular lines, and distributions of graphite and silicate 
grains. Decompositions of the SEDs of the template starburst galaxies NGC7714 and NGC2623 and the 
template AGNs PG 0804H-761 and Mrk463 provide baseline properties to aid in quantifying the strength of star- 
formation and accretion in the composite systems NGC 6240 and Mrk 1014. We find that obscured radiation 
from stars is capable of powering the total dust emission from NGC 6240. The presence of a small quantity of 
1260K dust in this source suggests a ^2% AGN contribution, although we cannot rule out a larger contribution 
from a deeply embedded AGN visible only in X-rays. The decomposition of Mrk 1014 is consistent with ^65% 
of its power emerging from an AGN and ~35% from star-formation. We suggest that many of the variations 
in our template starburst SEDs may be explained in terms of the different mean optical depths through the 
clouds of dust surrounding the young stars within each galaxy. Prompted by the divergent far-IR properties 
of our template AGNs, we suggest that variations in the relative orientation of their AGN accretion disks with 
respect to the disks of the galaxies hosting them may result in different amounts of AGN-heated cold dust 
emission emerging from their host galaxies. We estimate that 30-50% of the far-IR and PAH emission from 
Mrk 1014 may originate from such AGN-heated material in its host galaxy disk. In a subsequent paper, this 
decomposition method will be applied to a large sample of dusty galaxy SEDs to further test the validity of 
these suggestions. 

Subject headings: galaxies: Seyfert — galaxies: starburst — infrared: galaxies — methods: numerical 
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1. INTRODUCTION 

Optical and near-IR imaging has revealed that most Lu- 
minous and Ultra-Luminous Infrared Galaxies (LIRGs and 
ULIRGs — see Sanders & Mirabel 1996), as well as many 
lower luminosity dusty galaxies (e.g. starbursts and AGNs), 
are interacting systems which have recently merged or are 
in the process of merging (Armus et al. 1987; Sanders et al. 
1988; Murphy et al. 1996). Large quantities of gas and dust 
driven into the nascent nuclei of these systems catalyze mas- 
sive bursts of star formation and/or fuel accretion onto super- 
massive black holes (Mihos & Hernquist 1996; Struck 1999). 
Ultraviolet and optical emission from many of these dusty 
galaxies is significantly obscured by the same gas and dust 
providing their infrared power — a fact which greatly compli- 
cates the process of diagnosing the star-formation and accre- 
tion contributions to their bolometric luminosities. LIRGs and 
ULIRGs are believed to be the largest contributors to the far- 
IR background and star-formation energy density at z « 1-3 
(Blain et al. 2002; Elbaz & Cesarsky 2003), and local ULIRGs 
are frequently used as analogues for the luminous sub-mm 
galaxy population. The development of a proper understand- 
ing of the histories of star-formation and galaxy evolution in 
the universe therefore requires the construction of robust and 
reliable methods for use in diagnosing the dominant luminos- 
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ity source in obscured dusty galactic nuclei. 

Great progress towards this goal was made using data from 
the Infrared Space Observatory (ISO). Diagnostic diagrams 
based upon high-ionization mid-IR lines (e.g. [O IV] and 
[Ne V]), the slope of the mid-IR continuum, as well as the 
strengths of the 6.2 and 7.7 fim PAH emission features were 
constructed to quantify the relative contributions from ac- 
cretion and starburst activity in many sources (e.g. Genzel 
et al. 1998; Lutz et al. 1998; Laurent et al. 2000; Tran et al. 
2001). These observations provided evidence that ULIRGs 
are powered primarily by star-formation, with the fraction of 
AGN activity increasing with bolometric luminosity. How- 
ever, due to the limited wavelength coverage and sensitivity 
of the ISO spectrometers, this analysis was limited to a rela- 
tively small number of bright nearby sources. These observa- 
tional challenges have largely been abated with the improved 
wavelength baseline and sensitivity of the Infrared Spectro- 
graph^ (IRS; Houck et al. 2004) on the Spitzer Space Tele- 
scope (Werner et al. 2004). Large samples of starburst galax- 
ies (Brandl et al. 2006), AGNs (Hao et al. 2005a; Weedman 
et al. 2005; Schweitzer et al. 2006), local ULIRGs (Armus 
et al. 2007), and their high-redshift counterparts (e.g. Houck 
et al. 2005; Yan et al. 2005; Lutz et al. 2005) have been ob- 
served with the IRS. As with the previous generation of stud- 
ies with ISO, diagnostic diagrams have been constructed using 
the properties of high-ionization lines (with much improved 
sensitivity to place strict limits on [Ne V] emission), mid-IR 
spectral slopes (with a much broader mid-IR wavelength base- 

* The IRS was a collaborative venture between Cornell University and Ball 
Aerospace Corporation funded by NASA through the Jet Propulsion Labora- 
tory and the Ames Research Center 
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line), and PAH feature emission (now detected to z > 2). 

In addition to these observational diagnostics, many tools 
have been developed to model the SEDs of emission from 
dusty galaxies. Models containing a variety of grains-sizes 
and chemical compositions have been used to calculate the 
spectrum of emission from astronomical dust embedded in 
single (e.g. Desert et al. 1990; Li & Draine 2001; Draine et al. 
2007) and multiple (e.g. Dale et al. 2001; Li & Draine 2002) 
incident radiation fields. Additionally, numerous radiative- 
transfer models have been pubUshed, including several to cal- 
culate the emission from an AGN and its obscuring torus (e.g. 
Nenkova et al. 2002; Dullemond & van Bemmel 2005). Fi- 
nally, methods intended to decompose the SEDs of galaxies 
into a number of components have also been developed. For 
example, Klaas et al. (2001) present such a method and use 
it to decompose the far-IR SEDs of a sample of ULIRGs into 
contributions from several dust-modified blackbodies at dif- 
ferent characteristic temperatures. 

We present in this work a new method used to decompose 
the SEDs of galaxies into emission from populations of stars, 
an AGN accretion disk, PAH molecules, atomic and molec- 
ular lines, and distributions of thermally heated graphite and 
silicate grains at different fitted characteristic temperatures. 
Given the many observational and theoretical techniques to 
study dusty galactic nuclei described above, it is fair to ques- 
tion our decision to introduce yet another. We are motivated 
to do so for two principal reasons. First, while the radiative- 
transfer methods incorporating realistic dust models described 
above provide the most physically complete method of cal- 
culating the SEDs of dusty galaxies, these methods are also 
the most dependent on the assumed geometry of the source. 
Unfortunately, detailed spatial information to constrain these 
assumptions is available for only a handful of nearby galax- 
ies, so the parameter space in which to search for best-fitting 
models in agreement with observations is both sizable and de- 
generate. Second, the methods for calculating the emission 
from dusty galaxies described above are primarily intended to 
model starburst galaxies or AGNs. None of them, however, 
is particularly well-suited to calculating the emission from 
sources containing a starburst and an AGN, each of which 
may be extinguished by different levels of obscuration (see 
§5.4 for an expanded discussion). 

Our goal for this paper is therefore to convince the reader 
that our decomposition method is well-suited to analyzing the 
SEDs of composite sources, and is motivated by the desire to 
address: (1) What energy source is primarily responsible for 
powering the emission from dusty galaxies? (2) How homo- 
geneous are the SEDs of dusty galaxies which have similar 
AGN and starburst fractions? and (3) What geometries can 
explain any inhomogeneities in the SEDs of similarly pow- 
ered dusty galaxies? To meet this goal we apply our decom- 
position method to a small sample of well-studied galaxies, 
including the relatively unextinguished starburst NGC7714, 
the obscured starburst NGC2623, the quasar PG0804H-761, 
and the Seyfert 2 Mrk463. These fits provide templates by 
which to understand the properties of pure starburst galax- 
ies and AGNs — the constituents of composite sources such as 
many ULIRGs. Additionally, as examples of such compos- 
ite sources, we present fits to the LIRG NGC 6240 (believed 
to be powered primarily by star-formation) and the ULIRG 
Mrk 1014 (believed to be powered primarily by accretion). 

We begin in §2 with an overview of the spectral decomposi- 
tion method. In §3 we detail the method used to calculate the 
spectrum of emission from the various source components. In 



§4 and §5 we describe the adopted dust model and the method 
used to calculate the thermal emission from dust. In §6 and §7 
we describe PAH feature and atomic and molecular line emis- 
sion, respectively. In §8 we describe the sample of galaxies 
analyzed in this work, in §9 we present their decompositions, 
and in §10 we present detailed analysis of the results. The 
breadth of our conclusions are clearly limited by the size of 
our sample. In a subsequent paper our method will be applied 
to a much larger sample of starburst galaxies, AGNs, LIRGs, 
and ULIRGs obtained as part of the Spitzer GTO program. 

2. SPECTRAL DECOMPOSITION METHOD^ 
2.1. Decomposition Components 
The principal assumption in our decomposition method is 
that the ultraviolet to millimeter wavelength (10 A- 1000 ^m) 
SED of a dusty galaxy may be completely described in terms 
of emission from evolved stars in the disk of the host galaxy 
(ISRF^), populations of young stars in starbursts (SB), an 
AGN accretion disk (AGN), PAH molecules (PAHs), atomic 
and molecular lines (lines), and several distributions of ther- 
mally heated graphite and silicate grains (hot, warm, cool, 
and cold). In this paper the ISRF, SB, and AGN components 
are collectively referred to as 'source' emission components 
since they represent the primary illuminating sources heating 
dust within a galaxy. The hot, warm, cool, and cold com- 
ponents are collectively referred to as 'dust' emission com- 
ponents. Assuming these components form a complete basis 
describing the emission from dusty galaxies, the observed flux 
density may be expressed as' 

L source L fj^^^ 

where ^source and Ljust are the total apparent luminosities of 
the source and dust components derived from the fit, Dl is the 
luminosity distance to the galaxy, a, = L,/Lsource and ctj = 
Lj/Ldast are the fitted contributions to the total source and 
dust luminosities from each component, ihe fv= fv/ J fvdu 
are the normalized flux densities of each source or dust com- 
ponent, and fl^^' and f^}""' are the fitted PAH and emis- 
sion line flux densities. We include a maximum of four dust 
components which are typically heated to characteristic tem- 
peratures of fhot « 1400 K, fwarm « 200 K, fcool ~ 80 K, and 
Tcoid ~ 35 K. The actual temperature of each dust component 
is determined from the fitted magnitude of its illuminating ra- 
diation field energy density. 

Figure 1 presents an illustration of one possible geometric 
structure of a dusty galaxy within which to interpret the roles 
played by each decomposition component in equation (1). 
Panel (a) depicts the optically-thin (at mid-IR wavelengths) 
disk of the host galaxy from which evolved photospheric 
(ISRF) and cold dust emission emerges. Panel (b) shows an 
enlarged view of the central region of the galaxy consisting 
of an AGN and a starburst surrounded by screens of cool and 
cold dust. Panel (c) depicts the composite nucleus of the dusty 
galaxy, while the details of each nuclear component are shown 
in panel (d). The starburst is composed of SB emission from 

^ See http://isc.astro.cornell.edu/~jam258/cafe for more information about 
the method and the Continuum And Feature Extraction (CAFE) software. 

^ All names of decomposition components are printed in italic type 
throughout this work. 

' Note that fx units may be substituted for fi, units throughout this work 
and that variables accented with tildes are free parameters in the fitting. 
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Fig. 1. — Example geometry of a dusty galaxy illustrating the potential 
roles played by each decomposition component. Panel (a) depicts the host 
galaxy disk containing cold dust and knots of ISRF emission from evolved 
stars. Panel (b) shows an enlarged view of the core of the galaxy containing 
a nuclear region embedded within obscuring screens of cool and cold dust. 
Panels (c) and (d) depict the composite nucleus. The starburst contains SB 
emission from young stars, warm dust in the surrounding shells, and possibly 
cool and/or cold dust in spherical screens around each individual star-forming 
cloud. The AGN contains an AGN accretion disk, hot and warm dust in the 
obscuring torus, and warm and cool dust in narrow-hne-region (NLR) clouds. 

young stars, warm dust in the shells surrounding the young 
stars, and possibly cool and/or cold dust in spherical screens 
around each individual star-forming cloud. The AGN consists 
of emission from the AGN accretion disk, hot and warm dust 
on the illuminated and shaded sides of the dust clouds creat- 
ing the obscuring torus (or dust in illuminated and shielded re- 
gions of a 'classic' torus with a smooth dust distribution — see 
Fig. 6), and warm and cool dust in narrow-Une-region clouds 
(again from dust in illuminated and shaded regions). We em- 
phasize that our method does not require the definition of any 
specific geometry (i.e. we derive the contributions from each 
emission component without the need to define its proximity 
to any other component), and the actual geometry of a dusty 
galaxy may look quite different from the one presented. The 
illustration in Figure 1 should therefore not be taken as a lit- 
eral description of the geometry of a particular galaxy (since 
our method requires none), but instead is intended to provide 
a structure within which to visualize the multiple roles poten- 
tially played by each decomposition component. 

We note that star-formation in the disk of the host galaxy 
and in starburst regions are both explicitly accounted for via 
the ISRF and SB components. Furthermore, while Figure 1 
depicts the specific case of a nuclear starburst, the actual loca- 
tion of the star-formation (e.g. nuclear versus off-center) has 
no effect on the decomposition (assuming that any off-center 
emission is contained within the spectroscopic slits and photo- 



metric apertures used to construct the SEDs). For example, in 
the decomposition of an off-center starburst, the cool and cold 
decomposition components (or portions thereof) could be in- 
terpreted as belonging to off-center obscuring clouds. Note, 
however, that if a galaxy contains both a highly-enshrouded 
nuclear starburst and a less-enshrouded off-center starburst, 
our decomposition method cannot explicitly fit both compo- 
nents simultaneously. Instead, a mean optical depth some- 
where between the optical depths of the two individual com- 
ponents is obtained. Thus, care must be taken when interpret- 
ing the optical depths inferred from a globally integrated SED 
(this admonishment is not unique to our analysis, but is a gen- 
eral statement about the nature of spatially integrated data). 

2.2. Decomposition Algorithm 

Once assembled, the ultraviolet to millimeter wavelength 
SEDs of the galaxies in our sample are fit to the model of the 
total flux density from equation (1). The best-fit parameters 
are determined using a Levenberg-Marquardt least-squares 
routine to minimize 

= E^* [ft'''\h)-fT\h)] ' , (2) 

k 

where the wu are the weights applied at each wavelength 
(see §2.3), f^^ is the observed flux density, and is the 
modeled flux density from equation (1). For photometric data, 
we integrate the modeled flux density over the filter transmis- 
sion curve (if available) to properly calculate the broadband 
flux density. This is important for photometric bands contain- 
ing emission and absorption features (e.g. L-band photome- 
try containing the 3.3 /im PAH, 3.1 /xm water-ice, and 3.4^m 
HAC features — see §4.8 and §6), as well as those bands cov- 
ering spectral regions that change rapidly with wavelength 
(e.g. J-band photometry sampUng extinguished photospheric 
emission— see the near-IR fit to NGC 6240 in §9.1). 

Experience has taught us that the greater the number of 
free parameters in a fitting model, the less likely the result- 
ing 'best-fit' will converge to the desired global minimum in 
the function. Correspondingly, it becomes increasingly 
likely that the fit wiU converge to a local minimum, often with 
very unsatisfactory results. This is a practical matter and not 
a mathematical one since a fit with any number of parameters 
will reach the global minimum given sufficiently well-chosen 
initial conditions. In practice, however, it is desirable to have 
a more robust fitting method such that any reasonably chosen 
initial conditions will converge to the global minimum. To 
meet this goal, we have developed an algorithm which sep- 
arates the processes of fitting the continuum — defined as the 
sum of the source and dust components — and the individual 
PAH features and atomic and molecular emission lines. Such 
a division into two separate least-squares processes greatly re- 
duces the total number of free parameters in each fit. Param- 
eter errors are calculated after the first step of the fitting pro- 
cedure and propagated through to the following fit (see §2.3), 
so that the net accumulated uncertainty is reliably calculated. 

After an SED has been fit, we perform a series of tests to 
ensure that all components are well-constrained (i.e. that the 
parameter errors are smaller than their values). First, if the 
ISRF luminosity is < 1 % the sum of the luminosities of all 
source components, then the fit is repeated without the ISRF 
component. Similarly, if any of the observed dust component 
luminosities are < 1% the sum of the apparent luminosities of 
all dust components, the fit is repeated without those compo- 
nents. AH three source components and the hot dust compo- 
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nent emit strongly at ultraviolet to near-IR wavelengths. For 
some galaxies these overlapping emission components are de- 
generate and result in unconstrained parameters. Our second 
check is therefore to test if either the SB or AGN components 
are unconstrained (i.e. if their luminosity uncertainties exceed 
their luminosities). If the SB component is unconstrained, 
we fix the SB-to-PAHs luminosity ratio (i.e. the ratio of the 
unextinguished SB luminosity to the PAHs luminosity — see 
§9.2) to the value derived from fits to the template starbursts 
NGC7714 and NGC2623, and restart the fit. Similarly, if 
the AGN component is unconstrained, we fix the AGN-to- 
hot dust luminosity ratio (i.e. the ratio of the unextinguished 
AGN luminosity to the unextinguished hot dust luminosity — 
see §9.2) to the value derived from fits to the template AGNs 
PG 0804+761 and Mrk463, and restart the fit. As our final 
test, if any of the fitted optical depths are unconstrained, we 
fix the most unconstrained to zero and restart the fit. This se- 
quence provides a well-defined method by which to sequen- 
tially impose constraints until a successful fit is obtained. 

The steps in the fitting method are therefore: (1) Input the 
observed SED and the initial estimates of the PAHs and lines 
components (see §6 and §7); (2) Subtract the initial PAHs and 
lines components from the observed SED to create an ob- 
served continuum emission spectrum; (3) Fit this spectrum 
to obtain the model continuum; (4) Subtract the model con- 
tinuum from the observed SED to create an observed PAH 
feature and line emission spectrum; (5) Fit the individual fea- 
tures in this spectrum; (6) Perform the sequence of tests to 
ensure that the fit is well-constrained; and (7) If necessary, 
repeat the fitting process with the constrained components. 

2.3. Data Weighting 

The success or failure of a fitting algorithm and the reliabil- 
ity of the results is largely dependent upon the chosen method 
of weighting the data in equation (2). We adopt the function 



A, 



(3) 



where (j{Xk) is the total 1-cr uncertainty in the flux density at 
wavelength A^, and is a term which compensates for non- 
uniform wavelength sampling (i.e. ensuring that a region of 
the SED does not dominate the fit simply because it is sam- 
pled more finely — see Appendix A). The total uncertainty in 
any observed data is given by the quadratic sum of the statis- 
tical and calibration (i.e. systematic) uncertainties. For IRS 
spectroscopic data, we use the difference of two spectra ob- 
tained at different nod positions (as produced during a typical 
IRS observation) to calculate the statistical uncertainties (see 
Appendix B for details). Carefully extracted low-resolution 
IRS spectra may still exhibit some residual calibration uncer- 
tainties (e.g. fringes in the long-low modules), which we esti- 
mate to be ^2% of the flux density at any given wavelength. 
For photometric data, the total uncertainty is taken to be equal 
to the quadratic sum of the uncertainties provided in the litera- 
ture and an absolute calibration uncertainty (e.g. from scaling 
to the IRS apertures) taken to be 5%. The actual systematic 
uncertainties in the photometric data are difficult to estimate, 
although we note that the adopted uncertainties yield reason- 
able reduced values, so are likely good estimates. 

As described in §2.2, each decomposition consists of two 
separate least-squares fits. In the first (step 3 above), the flux 
density to fit is obtained by subtracting the PAHs and lines 
components from the observed SED, i.e. f^, = /^^'^-/^'^^^ — 



iSRf 

AGJV 



SB (2 Myr) 
SB (10 Myr) ffl 
SB (100 Myr) 




\ \ 
\ \ 
\ 

\ 



\ \ 



10"' 10"- 0.1 1 10 10" 

X (nm) 

Fig. 2. — Normalized lOA-lOOO/^m energy densities of the principal lu- 
minosity sources in dusty galaxies including: emission from starbursts 2, 10, 
and 100 Myr after an instantaneous burst (SB); evolved disk stars (ISRF); and 
an AGN accretion disk {AGN). 

yhnes \Yhen a component is subtracted from observed data to 
create an SED to be fitted, the uncertainty of the subtracted 
component must be propagated through the analysis. Thus, 
the total uncertainty input into the least-squares routine for 
this first fit is = CTj^ta + '''pahs + "'unes ■ Here, (Tdata is the un- 
certainty in the observed flux density, while ctpahs and CTunes 
are the uncertainties in the estimated PAHs and lines compo- 
nents, respectively. Since the estimated strengths of the PAHs 
and lines components are derived from the observed spec- 
trum, (TpAH,s and (Tiincs are dominated by the IRS uncertainties. 
We therefore adopt ctpahs = cunes = 0, since the uncertainties 
in the IRS data on which they depend are already included 
through (Tdata- In the second fit (step 5 above), the flux density 
to fit is obtained by subtracting the fitted continuum from the 
observed SED, i.e. = .Z^"'" -/^™', so that the total uncer- 
tainty input into the least-squares routine is = o'data + '^'conf 
Here, (icont is the formal uncertainty in the continuum calcu- 
lated using the full covariance matrix from the first fit. 



3. SOURCE EMISSION 

3.1. ISRF Component Emission 

We adopt the model of the average interstellar radiation 
field presented in Mezger et al. (1982) and Li & Draine (2001) 
to represent the SED of emission from stars in the disk of a 
galaxy. In this model, the mean energy density in the solar 
neighborhood, cu^^^^ = fl^^^'® = 4TrJlf^^'^, is given by 

A<912A, 
912A< A< 1100 A, 
1100A< A< 1340A, (4) 
1340A < A < 2460 A, 
En^nB.iT,) A > 2460 A, 

where is the Planck function per unit frequency and the 
W„ = [0.025, 0.25, 1] are weighting factors for each blackbody 
component of temperature T„ = [7500 K, 4000 K, 3000 K]. 
This model spectrum is shown in Figure 2, normalized to have 
unit integrated energy density. Emission from the ISRF de- 
composition component is assumed to emerge nearly unob- 
scured from a galaxy disk, so that the line-of-sight extinction 
to the emitting stars is small in the infrared (A > 1 /im). The 
flux density per unit frequency interval of emission from the 
ISRF component for use in equation (1) is therefore 
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Fig. 3. — Example SEDs of emission from a starburst (SB) lOand 100 Myr 
after an instantaneous burst, evolved disk stars (ISRF), and an AGN accretion 
disk (AGN). Also shown are the SEDs of the SB components extinguished 
by obscuring screens with the indicated optical depths, and the SEDs of the 
AGW component behind an Av ~ 25 screen for the indicated covering factors 
(dash-dotted lines). Note that emission from the 2 Myr starburst population 
declines steeply at A > 3000 A and that the ISRF component is very similar 
to a 3500 K blackbody for A > 1 fim {dotted lines — see text). 



An example ISRF component SED, fl^^^, is shown in Fig- 
ure 3. Displayed alongside the ISRF component is the emis- 
sion from a 3500 K blackbody (characteristic of evolved stars) 
scaled to match at 1 .5 /im. The two curves are nearly identical 
at A > 1 /zm, so that fitting near-IR data with the ISRF com- 
ponent is very similar to fitting it with a 3500 K blackbody. 

3.2. SB Component Emission 

In contrast to the modest ^IMQyr"' star-formation rate 
characteristic of emission from the ISRF component, the rates 
in starbursts are often much higher (e.g. > 10-100 MQyr"'), 
resulting in significantly different stellar populations. In or- 
der to properly characterize the spectral properties and bolo- 
metric luminosities of these populations, we use Starburst99 
(Leitherer 1990) to generate SEDs of instantaneous bursts at 
different epochs. We assume a Kroupa IMF having a M'^ 
power-law in the 0. 1-0.5 M0 range and a M'-^ power-law 
in the O.5-1OOM0 range. We further assume Padova AGB 
tracks with solar metallicity. Populations of stars in starbursts 
evolve dramatically during the first 10 Myr after a burst. At 
4 Myr the most massive stars still exist on the main sequence, 
but by 10 Myr all of the photoionizing stars have evolved into 
supergiants. As shown in Figure 2, the emission from a star- 
burst 2 Myr after a burst emerges almost exclusively at ul- 
traviolet wavelengths. Since we have very little data in this 
wavelength range for the sources in our sample, we are unable 
to explicitly constrain the presence or absence of such a pop- 
ulation. In contrast to this, the 10 and 100 Myr populations 
emit a large fraction of their luminosity at optical and longer 
wavelengths, where our SEDs are better constrained. There- 
fore, in order to adequately sample the evolutionary stages in 
a starburst that we are sensitive to, and to smoothly transi- 
tion into the evolved population modeled with the ISRF com- 
ponent, we create a composite starburst spectrum containing 
equal contributions from the 10 and 100 Myr populations- 



i.e. M 



SB . 



(m^*^'" + m^'^'"")/2. We note, however, that the omis- 
sion of an ultraviolet-luminous young stellar population from 
may result in an underprediction of the total SB luminos- 
ity for some sources, perhaps by up to ~l/3 if the 2 Myr 
population contributes at the same level as each of the older 



populations. 

Unlike emission from the ISRF component, we assume that 
SB component emission emerges from regions that may be ob- 
scured by screens of dust out to infrared wavelengths. If the 
obscuration to the SB component is dominated by dust in the 
cocoons of material surrounding each star-forming knot, and 
not by extinction between the various knots, then the majority 
of the obscuration occurs via a screened geometry. Leven- 
son et al. (2007) use radiative transfer models to show that 
such a screened geometry is actually required to reproduce 
the deep 9.7 /im silicate features observed in many ULIRGs, 
suggesting that the assumption of a screened geometry is rea- 
sonable for highly obscured sources. At lower optical depths 
(i.e. T < 1), mixed and screened geometries both give rise 
to extinctions scaling approximately as ///o ~ 1 - r, so that 
the choice between the form of extinction is less critical. We 
therefore model the flux density per unit frequency interval of 
emission from the SB component as 

/fcxM^e-^-, (6) 

where the optical depth through the obscuring screen to the 
SB component is given by 

^ / \\ _ ^SB^^5abs(A)^ 

TSB(A) = Ty —. (/) 

I]ab.s(5500A) 

Here, fy^ is the screen optical depth to the SB component 

at 5500 A, and Eaks = Sext- Sjcat is the total dust absorption 
opacity obtained by subtracting the scattering opacity from 
the total extinction opacity (see §4.7). We use the absorption 
and not the extinction opacity in equation (7) since scattering 
tends to cancel itself out in sufficiently spherical geometries, 
such as depicted in Figure 1 (i.e. light scattered out of the 
beam along the direct path to the observer is offset by light 
scattered into the beam along another path). Furthermore, 
the extinction and absorption opacities of dust in our model 
have similar slopes between 1000 A and 3 /im, where the SB 
component emits strongly. Since the derived optical depth de- 
pends only on the opacity slope and not its absolute value, the 
omission of scattering has very little effect. We note, how- 
ever, that this omission may introduce some uncertainties at 
A < 3 /im, where the scattering opacity is largest. Example 10 
and 100 Myr SEDs, Z^'^"' and Z^'^'"", used to construct the 
composite SB model, are shown in Figure 3. Also shown 
are example SEDs obscured by screens of dust with f^^ = 1, 
2, and 4. 

3.3. AGN Component Emission 

The spectral properties of radiation emerging from an AGN 
accretion disk depend sensitively upon the assumed geometry 
of the surrounding obscuring structure — i.e. the putative torus 
of AGN unification models (Antonucci 1993). Important ge- 
ometrical factors include the height of the inner edge of the 
torus, its proximity to the nucleus, the clumpy versus smooth 
structure of the obscuring medium, and the orientation of the 
torus with respect to the observer. Radiative-transfer calcula- 
tions incorporating these geometric properties have been used 
to construct models of the emission from an AGN accretion 
disk and its obscuring torus (e.g. Pier & Krolik 1992; Laor & 
Draine 1993; Granato & Danese 1994; Nenkova et al. 2002; 
Dullemond & van Bemmel 2005). Since our intent in this 
work is to characterize the general properties of the emergent 
SEDs of large samples of galaxies, and not to discern the de- 
tailed properties of the obscuring torus in individual sources, 
such modeling is unnecessary for our purposes. 
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Instead, we adopt an empirical model of the emission from 
an AGN accretion disk and assume that this SED may be ex- 
tinguished by a foreground screen of dust. We use the accre- 
tion disk model from Schartmaim et al. (2005) in which the 
emergent energy density, shown in Fig. 2, is modeled by 





A3 

1. 5^(1000 K) 



A < lOA, 
lOA < A < 500 A, 
500A< A< 1216A, 
1216A< A< 10/xm, 
A > lOyum. 



(8) 



The flux density per unit frequency interval of emission from 
the AGN component, including an obscuring screen with a 
non-uniform covering factor, is modeled as 

oc [(1 - eAGN) + eAGNe-"^™] , (9) 

where < cagn < 1 is the fraction of the accretion disk ob- 
scured by a screen of optical depth tagn, with 

Sext(A) 



tagn(A) = fy 



-AGN 



Sext(5500A) 



(10) 



Here, Ty™ is the screen optical depth to the AGN component 
at 5500 A. An example unobscured AGN component SED, 
Z^'^^, is shown in Figure 3. Also shown are the SEDs of AGN 
component emission obscured by a fy'^^ = 25 screen, for sev- 
eral values of the screen covering factor eACN- Note that sig- 
nificant ultraviolet and optical emission may still emerge from 
a highly extinguished accretion disk if the covering factor is 
less than unity. 

In the bottom panel of Figure 1 , the AGN accretion disk is 
depicted as being surrounded by an obscuring structure com- 
posed of discrete clouds and/or a smooth dusty component. 
The covering factor, cagn, may therefore be interpreted as the 
fraction of the accretion disk covered by this obscuring ma- 
terial (from the point-of-view of the observer), which has an 
approximate visual band optical depth equal to f^*^'^. Note 
that the same qualifications regarding the relative importance 
of scattering and the subsequent uncertainties in the emergent 
flux density for A < 3 /xm that were made for the SB compo- 
nent apply here as well. Since the clumpy torus is by def- 
inition not spherically symmetric, we use the full extinction 
opacity (including scattering) in equation (10). As described 
for the SB component, this choice has very little effect on the 
results since the extinction and absorption opacities have sim- 
ilar slopes from the ultraviolet to the near-IR, where the AGN 
component emits strongly. 

4. ASTRONOMICAL DUST MODEL 

4.1. Grain Properties 

We adopt the grain-size distribution function from Wein- 
gartner & Draine (2001, hereafter WDOl), and the optical 
properties of graphitic carbon and smoothed astronomical sil- 
icate from Draine & Lee (1984) and Laor & Draine (1993), 
with modifications from Weingartner & Draine (2001) and Li 
& Draine (2001). The WDOl distribution function has been 
modified to include a small grain- size exponential cutoff to 
model the deficit of small grains resulting from sublimation 
in intense heating conditions (see §4.2). The modified distri- 
bution function for graphite and silicate grains takes the form 



1 dni 
hh da 



1 dni 
hh da 



1-exp 



WDOl 



13-' 



(11) 



where [«jj'(i«,/t/a]wDOi is the WDOl distribution function for 
grains with radii satisfying 3.5A < a < 1 /im, = a^(U) is 
the minimum surviving grain-size in a distribution embedded 
in a radiation field with an energy density of magnitude U (see 
§4.2 and eq. [12]), and /?_ determines the exponential cutoff 
rate. In the limit of oo, the exponential provides a sharp 
cutoff at a = a-. We choose /3_ = 3, the same value adopted for 
the large grain-size exponential cutoff in WDOl. Absorption 
and extinction efficiencies, Q'^^^^ and gexf h^^^ hem tabulated 
for graphite and silicate spheres as described in Laor & Draine 
(1993) for wavelengths lOA < A < 1000/um. Absorption 
cross sections are calculated using C^^(a,A) = 7ra^(2^s(a, A) 
and C^b/fl, A) = Tra^Qf^^(a, A), with similar expressions for the 
extinction cross sections Cf^t (a, A) and QJ;\(a, A). 

4.2. Grain Sublimation 

Graphite and silicate grains sublimate when heated to 
temperatures above 1750 and 1400 K, respectively (Laor & 
Draine 1993). The minimum surviving grain-size in a dis- 
tribution, a-, is determined by the size of the smallest grain 
which does not sublimate in the radiation field in which the 
distribution is embedded. As the magnitude of the radiation 
field energy density increases, the grain distribution becomes 
increasingly depleted of small grains. A stronger radiation 
field therefore changes the SED of emission from a distri- 
bution of grains both by increasing the temperature of the 
grains, and by weighing the grain-size distribution function 
to larger grains. Additionally, due to their lower sublimation 
temperatures, silicate grains become depleted before graphite 
grains. SubUmation is one of many physical processes which 
can set the minimum grain-size within a distribution. Shock 
waves (Jones et al. 1994; Draine 1995), sputtering (Draine 
& Salpeter 1979), and grain-grain collisions (Tielens et al. 
1994) can also influence the minimum grain-size. However, 
these other effects are more difficult to model than subUma- 
tion, since they require a dynamic model of the region where 
the grains are located. We do not consider them further. 

4.3. Thermal Heating 

The equilibrium temperature attained by a grain is a func- 
tion of its composition, size, and the spectrum and strength of 
the radiation field in which it is embedded. Our adopted dust 
model contains grains spanning several orders of magnitude 
in size, so the temperatures of individual grains within a dis- 
tribution embedded in a given radiation field vary greatly. For 
a graphite or silicate grain embedded in a radiation field with 
normalized energy density Hi, (see Fig. 2) which has a fitted 
magnitude U, energy conservation requires that 



Ci^,,(a,iy)Ucu^du= j C[^^(a,u)ATiB,[Tl^{a)]dv. (12) 



Here, the left and right sides represent the rates at which en- 
ergy is absorbed and emitted by a grain, respectively. Equa- 
tion (12) is therefore an implicit equation for the equihbrium 
temperatures, T^^ia), of the graphite and siUcate grains in a 
distribution. These temperatures are used in equation (15) to 
properly calculate the spectrum of emission from each de- 
composition component. We calculate the temperatures of 
grains in the hot dust component using the normahzed AGN 
energy density, m^™, and the temperatures of grains in the 
warm, cool, and cold dust components using the normalized 
ISRF energy density, u^^^^ (see Fig. 2). Note that although 
the shape of the iUuminating radiation field does affect the 
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temperatures of grains (since grains absorb and emit more ef- 
ficiently at certain wavelengths), similar grain temperatures 
may be obtained from two different radiation fields by suit- 
ably adjusting the magnitude of the energy density of each 
field. Thus, the specific choices about which illuminating ra- 
diation fields we use to calculate our dust equilibrium temper- 
atures do not have a significant impact on our results. Note 
also that once the shape of the normalized energy density is 
chosen, the temperature of a grain in that field will scale with 
the fitted magnitude of the energy density, U . 

4.4. Stochastic Heating 

The heat capacities of very small grains are sufficiently 
small that they undergo large temperature fluctuations upon 
absorption of individual photons, and reradiate most of the 
deposited energy on timescales much shorter than the typical 
time between photon strikes. To estimate the importance of 
this stochastic heating for our dust components, we calculate 
the threshold grain-size below which grains fall out of thermal 
equilibrium. For graphite grains, Draine & Li (2001) calcu- 
late this to occur for 



amin ~ lOOA f 



\Mmmp 



-0.2 



(13) 



where u = J Ui,dv is the integrated spectral energy density 
of the illuminating radiation field, and mmmp is the integrated 
spectral energy density of the local interstellar radiation field 
from Mezger et al. (1982). For our cold, cool, warm, and 
hot dust components, typical threshold grain-sizes from equa- 
tion (13) are Omin ~ 53, 14, 7.5, and 3 A. 

In the model of Li & Draine (2001), stochastically heated 
grains exhibit strong spectral feature emission (i.e. PAHs) but 
have very weak continua (due to the low continuum opaci- 
ties of small grains). Most dusty galaxies contain a significant 
amount of thermally heated f > 1 00 K dust which emits more 
continuum radiation than these stochastically heated grains 
(see also Laor & Draine 1993). Thus, the continuum radia- 
tion from stochastically heated grains in our cold component 
(i.e. grains with a < 53 A) is weak compared to the thermal 
emission from warmer dust. Additionally, since the threshold 
grain-size decreases as the magnitude of the illuminating radi- 
ation field energy density increases, stochastically heated con- 
tinuum emission is even less important for grains in our cool, 
warm, and hot components. Since the stochastic continuum 
emission of all components is weak compared to the emission 
from thermally heated grains, and given that the stochastically 
heated feature emission attributed to PAHs is included in our 
model through the PAHs component, we do not incorporate a 
detailed model of stochastic heating. 

4.5. Dust Emissivity 

The total emissivity per hydrogen nucleon of a distribution 
of grains (in units of erg s"' Hz"' sr"' H"') is obtained by 
summing over the emissivity of each grain species 



Ely — E^ + Ej^ 



(14) 



Here, the emissivities of graphite and silicate grains are 

Ei. 



1 dn 

^ —^Ci^„(a,X)BAT'^(a,f)]da, (15) 
50 A «H da 

which is an average over the individual emissivities of the 
grains in the distribution. Note that in contrast to models 



100 : 



Cold 



T(K) 

Fig. 4. — Range of equilibrium temperatures, Tcq, of tlie grains within 
a distribution as a function of the characteristic distribution temperature f. 
The dark shaded region shows the range of equilibrium temperatures within a 
given distribution (i.e. at a fixed value of T), while the solid line shows their 
luminosity-weighted characteristic temperature. Light shaded regions indi- 
cate typical characteristic temperature ranges for the decomposition compo- 
nents. The sublimation temperature of silicates is also labeled {dashed line). 

using a single effective grain-size, the dust emissivity in our 
adopted model is not directly proportional to the dust opacity 
(see §4.7), since the equilibrium grain temperatures in equa- 
tion (15) depend upon grain-size and composition. 

4.6. Characteristic Distribution Temperature 

As discussed in §4.3, the many grains within a distribution 
are brought to different equilibrium temperatures through in- 
teractions with the radiation field in which they are embedded. 
In our decomposition method this distribution of equilibrium 
temperatures is determined from the fitted magnitude of the 
illuminating radiation field energy density, U . A given dust 
component therefore has many temperatures associated with 
it, none of which uniquely characterize the properties of the 
distribution. This is clearly not a desirable situation, as we 
would like to associate a single characteristic temperature to 
each value of U, preferably a temperature which communi- 
cates the position of the spectral peak of the resulting emis- 
sion from the distribution. To meet this goal, we define the 
luminosity-weighted characteristic distribution temperature 



T = T(U) = 



I^Aj:,fL(^,U)T'^(a,U)da 



(16) 



I^A^JL(^,U)da 
where the sum is over the two grain compositions, and 

fi(a,U)=—^ f Ci^,,(a,iy)BAT:^(a,U)]diy (17) 
nu da J ^ 

is the grain-luminosity distribution function (which is es- 
sentially the spectrally integrated emissivity of an individual 
grain — see §4.5). Note that the integration starts at a = 50 A, 
since smaller grains are not in therrnal equilibrium. 

Using this weighting function, f is properly interpreted 
as the composition-averaged temperature of the luminosity- 
dominating grain-size. With this definition, the emission peak 
from a distribution of grains occurs near the wavelength ex- 
pected from the dust-modified Wien's law 



A, 



peak ■ 



(^) 



lO^m, 



(18) 



where the peak is defined to occur at the wavelength coincid- 
ing with the maximum value of E^y = Y,^i,sB,^(T) (i.e. the emis- 
sivity of a distribution in the limit in which the temperatures of 
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Fig. 5. — Water-ice and HAC opacity templates derived from L-band 
(Imanislii et al. 2006) and IRS (Spoon et al. 2004) spectra of ULIRGs (see 
text). Also shown in the left panel is the arbitrarily scaled IRAC 3.6 /^m 
transmission curve for sources at z = 0, 0.1, and 0.2 (from right to left). 

its constituent grains are independent of their size and compo- 
sition). Figure 4 shows the range of equilibrium temperatures, 
Teq, of the grains within a distribution for a given value of the 
characteristic distribution temperature, f . The dark shaded 
region spans the range of equilibrium temperatures in a dis- 
tribution, while the solid line shows the luminosity-weighted 
average over that distribution. Note that this average tempera- 
ture is weighted towards the warmer equilibrium temperatures 
since the grain-size distribution function is weighted towards 
smaller (i.e. warmer) grains. The light shaded regions indi- 
cate the typical ranges of the characteristic distribution tem- 
peratures for the dust components in our decomposition (i.e. 
as determined by the fitted values of I] for each component). 

4.7. Dust Opacity 

The total absorption opacity per hydrogen nucleon of a dis- 
tribution of grains (in units of cm^ H"' ) is obtained by sum- 
ming over the opacity of each grain species 



^abs(A)=EfrsW+Sabs(A). 



Here, the opacities of graphite and silicate grains are 



SL(A) = 



^ 1 dn; : 
50 A "H da 



(19) 



(20) 



which is an average over the individual opacities of the grains 
in the distribution. The total extinction opacity per hydrogen 
nucleon, Eext, is calculated in an analogous way with Cab.s ^ 
Cext. According to Li & Draine (2001, Table 6), Sabs oc 
over the range 20 [im < A < 700 /im and Eabs oc A" 
the range 700 ^m < A < 10''/^m to within ±10%. 



-1.68 



4.8. Water-ice and HAC Opacity 

Mid-IR spectra often exhibit 5.6-7.8 pm opacity produced 
from a combination of water-ice absorption near 6.1 pm (see 
e.g. Chiar et al. 2000; Gibb et al. 2000) and absorption from 
hydrogenated amorphous carbon (HAC) near 6.85 pvs\ and 
7.25 /im (Furton et al. 1999). Such water-ice and HAC ab- 
sorption is particularly prominent in ULIRG spectra (e.g. 
Spoon et al. 2002). We derive 5.6-7.8 pm water-ice and HAC 
opacity templates from the IRS spectrum of the heavily ob- 
scured ULIRG IRAS F00183-711 1 (Spoon et al. 2004). This 
high signal-to-noise spectrum exhibits strong water-ice and 
HAC absorption with very little R4H emission, thereby pro- 
viding a clean spectrum from which to extract profiles of 
the various opacity sources. We fit a smooth spline to the 
observed spectrum between continuum points on either side 



of the water-ice and HAC features, and use the ratio of the 
observed spectrum to the estimated continuum to obtain the 
opacity templates displayed in the right-panel of Figure 5. 

Additional opacity from 3.1 pm water-ice and 3.4/im HAC 
features may significant affect the broadband photometry of 
some sources (see Fig. 5). We therefore use the L-band spec- 
tra of ULIRGs presented in Imanishi et al. (2006) to de- 
rive opacity templates for these features to use in our de- 
compositions. The opacities of the 3.1 /xm water-ice and 
3.4/im HAC features are derived from spectra of the deeply 
obscured ULIRGs IRAS 00188-0856 and IRAS 08572H-3915, 
respectively. Note that we do not derive the 6.1 pm water- 
ice template from the IRS spectrum of IRAS 00188-0856 
since it is contaminated by 6.2 /im RA.H emission. As 
measured from their L-band and IRS spectra, the water-ice 
features of IRAS 00188-0856 have appai-ent T^^\lrf \ w 3.9 
while the HAC features of IRAS 08572H-3915 have apparent 



^HAC 



HAC , 
85 



: 3.3. Our full opacity templates, shown in Fig- 



ure 5, are constructed using these optical depth ratios to set 
the scalings between the 3. 1-3.4 pm and 5.6-7.8 pm features. 

4.9. Total Optical Depth 

The total optical depth along a line-of-sight is the sum of 
the optical depths from each opacity source. Thus, in our case 



t(A) = Tdust(A) + Tice(A) + THAc(A). 



(21) 



The optical depth through a column of dust having total ab- 
sorption or extinction (depending on the situation) opacity 
I](A) is given by 



Tdust(A)=A^HS(A) = f9^7 



S(A) 



E(9.7//m)' 



(22) 



where A^h and fg"^' are the hydrogen nucleon column density 
and 9.7 pm optical depth through the dust, respectively. The 
optical depth due to ice absorption is given by 



/ \ \ _ ~ -dust 
'icel^^ — '/ice '9 7 



Sice(A) 
Eice(6.1/im)' 



(23) 



where Eice is the ice opacity template from Figure 5, and 
Tjice = T^^\ I Tg determines the amount of ice opacity for a 
given dust opacity. Similarly, the optical depth due to HAC 
absorption is given by 



rHAc(A) = r/HACTg.T' 



Shac(A) 
SHAc(6.85//m)' 



(24) 



where Shac is the HAC opacity template from Figure 5, and 
77HAC = ■^6^85^ /■^gT' determines the amount of HAC opacity for 
a given dust opacity. The ratios of the apparent 6.1 pm water- 
ice and 6.85 pm HAC optical depths to the apparent 9.7 pm 
silicate optical depth are both ~0.1 for the various ULIRGs 
used to derive the opacity templates. We therefore restrict 
the fitted values of these ratios to satisfy < Tyice < 0.1 and 
< T/HAC < 0. 1 in our decomposition method. 

5. DUST EMISSION 

5.1. Optically -Thin Dust Emission 

Figure 6 presents an illustration of the optically-thin dust- 
shell approximation used to calculate the SEDs of emis- 
sion from the dust components in our decomposition that 
are heated directly by ultraviolet and optical photons from 
the SB and AGN source components (see also Fig. 1). The 
left-panel depicts a starburst containing a SB source compo- 
nent surrounded by a shell of dust. The right-panel shows 



Decomposing Dusty Galaxies. I. 



9 



Starburst 




AGN 

Clumpy -or- Smooth 

Torus Hot & Warm Torus 





NLR Clouds 
Warm & Cool 



Fig. 6. — Cross-sectional views illustrating the optically-thin dust-shell ap- 
proximation used to calculate the SEDs of dust components that are directly 
heated by ultraviolet and optical radiation from the SB and AGN source com- 
ponents (see text). The left-panel depicts a SB source component suiTounded 
by a shell of dust. The right-panel shows an AGN accretion disk and torus 
(made either of discrete clouds as on the left or a smooth dust distribution as 
on the right), and more distant narrow-line region clouds. 

an AGN comprised of an AGN accretion disk, a dusty torus 
(either smooth or made of discrete clouds), and more distant 
narrow-line region clouds. In this geometry, ultraviolet and 
optical photons emitted by the source components either es- 
cape freely or are absorbed within shells of dust on the illu- 
minated sides of clouds. By definition, these shells have an 
optical depth of order unity at the wavelength of peak energy 
absorption 

^ ^ _ j Y..^h^(\)ux\d\ ^^^^ 

Here, the absorption-weighted distribution function is ob- 
tained by integrating the energy absorbed by a single grain 
(i.e. the left-side of eq. [12]) over the distribution. As de- 
scribed in §4.3, we assume that grains in the hot compo- 
nent are heated by the AGN radiation field (Fig. 2), while 
grains in the other dust components are heated by the ISRF 
field. The peak absorption wavelengths from equation (25) 
for these fields are 1470 and 4550 A, respectively. At these 
wavelengths, the dust absorption opacity is approximately 20 
and 5 times greater than its value at 9.7 /^m, so that the 9.7 /im 
optical depths of the thin dust shells illuminated by the AGN 
and ISRF radiation fields are r^''^" w 0.05 and 0.2, respec- 
tively. 

All emission from these shells therefore emerges from re- 
gions which are optically-thin at all infrared wavelengths (i.e. 
A > 1 fjm). The emerging intensity from these shells, calcu- 
lated from the radiative transfer equation, is therefore 



-[Tshdl- 



Ut' = [l-e-^*"]5^. 



(26) 



Here, the source function S^, = fi^/Eabs, is assumed to be 
constant (i.e. the magnitude of the radiation field energy 
density — and thus the dust temperature distribution — is as- 
sumed to be uniform). As shown above, Tsheii < 1 for all in- 
frared wavelengths, so that equation (26) simplifies to 

/,«A^^^^"£,, (27) 

where A^^'"'" is the hydrogen nucleon column density through 
the dust shell (i.e. Tsheii = A^h ''"^abs)- 

5.2. Hot and Warm Component Emission 

Using equation (27), the flux densities per unit frequency 
interval of emission from the hot and warm components are 

fl°' cx £,(fhot)e-""°' and cx £,(fwa™)e"""™ , (28) 



B 0.10 




X (|lm) 

Fig. 7. — Example SEDs of emission from the hot, warm, cool, and cold 
dust components for typical characteristic temperatures of T = 1400, 200, 80, 
and 35 K (black lines). Also shown are the SEDs obscured by screens (gray 
lines) having rg 7 = 1,2, and 3 (from top to bottom). 

where 7i,ot = f{U^od and f„ai-m = r(J/„ai-m) ai-e the characteris- 
tic component temperatures determined from the fitted values 
of f/hot and {/warm (scc §4.6), and 



Thot ■ 



'9.7 



Sabs (A) 

Eabs(9.7/im) 



and r„ 



'97 



I Sabs (A) 

Sabs(9.7Mm) 



(29) 



are the optical depths through the screens of dust obscuring 
the hot and warm components (e.g. provided by the cool and 
cold components in the geometry of Fig. 1). Here, fg°' and 
f ™™ are the screen optical depths to the hot and warm com- 
ponents at 9.7 jivci. Note that we use the absorption and not the 
extinction opacity since scattering is negligible for A > 3 /im, 
where essentially all dust emission is radiated. Example hot 
and warm component SEDs, /^°' and Z*"™, are shown in Fig- 
ure 7 for the characteristic temperatures Thot = 1400 K and 
Twarm = 200 K. The fitted hot component characteristic tem- 
perature is constrained to satisfy 500 K < fhot < 1500K and 
the fitted warm component characteristic temperature typi- 
cally satisfies 150K < f„aim < 500 K. Also shown in Fig- 
ure 7 are example hot and warm component SEDs obscured 
by screens having tq 7 = 1,2, and 3. Note that departures from 
a smooth continuum near 9.7 and 18/im are caused by emis- 
sion and absorption from silicates. 

5.3. Cool and Cold Component Emission 

In the geometry of Figure 1, dust in the cool and cold com- 
ponents may serve as sources of opacity to obscure the emis- 
sion from other components. Given this geometry, and the 
observation that the mid-IR obscuration in dusty galaxies is 
often quite high (e.g. rg 7 w 3-5 for many of the local ULIRGs 
in Armus et al. 2006), the optically-thin approximation used 
to obtain equation (27) is not necessarily valid for the cool 
and cold components. However, as a result of their rela- 
tively cool characteristic temperatures and 40-100 /im SED 
peaks, the optical depth through these components at the 
wavelengths over which they emit strongly is a factor ^10- 
30 times smaller than at 9.7 /im. Therefore, despite poten- 
tially large 9.7 /im optical depths, both components are actu- 
ally optically-thin at the wavelengths over which the majority 
of their emission is radiated. Thus, the conditions required for 
the validity of equation (27) are in fact still satisfied so that the 
flux densities per unit frequency interval of emission from the 
cool and cold dust components are 



/»°' cx £,(7;ooi) and /™'" a El°'\T,,,i) 



■cold 



7 cold/ 



(30) 
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where Tcooi = T{Ucooi) and T^oid = T{Uco\d) are the character- 
istic component temperatures determined from the fitted val- 
ues of f/cooi and t/coid- We assume that emission from both 
the cool and cold components emerges unextinguished since 
the optical depths required to achieve significant extinctions 
at A > 50 /im where they emit strongly would imply implau- 
sibly high optical depths at shorter wavelengths. Example 
cool and cold component SEDs, and /^'^ are shown 
in Figure 7 for the characteristic temperatures fcooi = 80 K and 
Tcoid = 35 K. The fitted cool and cold component character- 
istic temperatures typically satisfy 50 K < fcooi < lOOK and 
25 K < fcoid < 40 K, respectively. Also shown in Figure 7 are 
the cool and cold SEDs obscured by screens having rg 7 = 1, 
2, and 3, demonstrating that the obscured cool and cold com- 
ponent SEDs are changed very little. 

5.4. Comparison with Other Methods 

The emergent flux density from dust embedded in a distri- 
bution of radiation fields with energy densities U is 

' dMd, 



jrdust 



dU 



-EAT(U)]dU. 



(31) 



Here, Mdust(t^) is the mass of dust heated to the character- 
istic temperature f(U) by the field with strength U. Our 
four-component dust model characterizing the emission from 
a galaxy is derived from equation (31) by taking 



dU 



(32) 



where the sum is over four discrete values of U, illuminating a 
mass M^uj,j of dust each. Clearly, equation (32) is an oversim- 
plification of the true distribution of mass in a galaxy. Dale 
et al. (2001) and Li & Draine (2002) both calculate the inte- 
grated emission from dust in galaxies by assuming a power- 
law distribution of dust mass over heating intensity 



dU 



(X U~ 



(33) 



where a determines the relative contributions from each value 
of the radiation field energy density. Given a sufficient num- 
ber of discrete temperature components, the model in equa- 
tion (32) can be used to approximate the emission from the 
power-law model to arbitrary precision (with the appropriate 
weighting of each component through the M'^^^^). The fact that 
both our multi-component model and these power-law mod- 
els produce excellent fits to the integrated spectra of galaxies 
implies that the required 'sufficient number' of discrete com- 
ponents postulated above must be approximately four. 

The real power of our multi-component approach stems 
from the fact that it is capable not only of approximating the 
emission from the mass distribution of star-forming galaxies 
in equation (33), but also the emission from mass distributions 
which are not characterized by a power-law (since we do not 
make the geometrical assumptions concerning the distribution 
of dust within a galaxy implied by eq. [33]). For example, 
the composite emission from a galaxy containing both a star- 
burst and an AGN may have a total mass distribution function 
dominated by a star-formation driven power-law at low energy 
densities superposed with an AGN driven power-law (having 
a different exponent) at high energy densities. Since the sum 
of two power-laws with different exponents is not necessar- 
ily itself a power-law, a model such as equation (33) cannot 



1.00 
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0.10 



0.01 




X (fim) 



Fig. 8. — PAH emission template (black line) derived from the spectrum 
of tlie mean starburst galaxy from Brandl et al. (2006). The individual com- 
plexes (gray lines) are modeled with Drude profiles (see Smith et al. 2007). 
Note the significant 'continuum' emission created by the addition of flux in 
the wings of the broad features. Also shown for reference are the four arbi- 
trarily scaled IRAC transmission curves. 

properly model this scenario. This, coupled with the fact that 
our method allows us to fit for different levels of obscuration 
towards each component (as expected for physically distinct 
AGN and starburst regions), drives our decision to utilize the 
multi-component discrete temperature approach to describe 
the composite SEDs of dusty galaxies. 

5.5. Dust Emission in Previous Publications 

The decompositions of IRAS 10214H-4724 in Teplitz et al. 
(2006), NGC6240 in Armus et al. (2006), and the BGS sam- 
ple of ULlRGs in Armus et al. (2007) were performed using 
an earlier version of the method described in this work which 
utilized a different model to calculate the emission from each 
dust component. In these previous publications, dust emis- 
sion is calculated from an optically thin shell of constant den- 
sity material surrounding the illuminating source. As in equa- 
tion (33), the dust in this shell is exposed to different radia- 
tion fields (depending on its distance from the central source), 
and is therefore brought to different equilibrium temperatures. 
In these previous works, we assume that the dust shells had 
a thickness rout/nn =10, which fixes the range of radiation 
fields and therefore dust temperatures, once the temperature 
at the minimum radius is defined (taken to be the sublimation 
temperature). Since each of these dust components contains 
a range of temperatures (from dust at different radii), only 
three components were needed to fit the sources instead of 
the four used in the current work. We choose to adopt the 
four-component method for this work since it improves the 
fits at long wavelengths and simplifies the interpretation. 

6. PAH EMISSION 

Observations with ISO and Spitzer have shown that the 
presence of PAHs (Leger & Puget 1984; Allamandola et al. 
1985) can be used as an indicator of star-formation (e.g. 
Forster Schreiber et al. 2004; Peeters et al. 2004). We note, 
however, that Calzetti et al. (2005) find that PAH emission 
does not correlate well with other star-formation tracers in 
M51. Additionally, Engelbracht et al. (2005) and Wu et al. 
(2006) find that PAH emission varies with metallicity, so 
may therefore not be a useful tracer of star-formation in low- 
metallicity systems. With these cautionary notes in mind, the 
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TABLE 1 
PAH Template Parameters 



TABLE 2 

Atomic and Molecular Line Parameters 



rCtiturc 
j 


^0 

(urn) 


7i 
(10-2) 


J 1^,0 

(Jy) 


r . / J .a,b 

(10-2) 


L-Ornplcx 

j 


I 


3.30 


1.2 


0.07 


0.6 


1 


2 


5.27 


3.4 


0.06 


1.0 


2 


3 


5.70 


3.5 


0.07 


1.1 


3 


4 


6 22 


n 
j.u 


77 


Q 7 


4 


5 


6.69 


7.0 


0.10 


2.7 


5 


5 


7.42 


12.6 


0.24 


39.3 




7 


7 60 


4 4 


89 








7 85 


5.3 


86 






9 


8 33 


5 


25 


3 9 


7 


10 


8 61 


3 9 


65 


7.6 




11 


11.23 


1.2 


1.29 


11.3 


9 


12 


11.33 


3.2 


1.04 






13 


11.99 


4.5 


0.32 


3.2 


10 


14 


12.62 


4.2 


0.85 


7.9 


11 


15 


12.69 


1.3 


0.19 






16 


13.48 


4.0 


0.22 


1.7 


12 


17 


14.19 


2.5 


0.18 


0.8 


13 


18 


16.45 


1.4 


0.43 


8.4 


14 


19 


17.04 


6.5 


0.64 






20 


17.38 


1.2 


0.29 






21 


17.87 


1.6 


0.26 






22 


18.92 


1.9 


0.31 


0.8 


15 



^ Fraction of total PAH luminosity emerging from eacii complex j. Lj = 
^.^^ L, is tile total luminosity of a complex of Drude protiles, wiiere L, = 
47rD^(7r/2)AJj7,/^^ Q is tiie apparent total luminosity of a single profile, 
with Dl the luminosity distance to the source.'^ As defined in Smith et al. 
(2007), the 7.7, 11.3, 12.7, and 17 /xm complexes are each composed of 
multiple Drude profiles. 

presence or absence of PAH emission nevertheless provides 
an important clue to help distinguish between star-formation 
and accretion dominated ULIRGs (e.g. Genzel et al. 1998), 
thereby necessitating the development of a method to care- 
fully decompose the contributions from PAHs to the total 
emission from dusty galaxies. 

This process is complicated by the fact that it is difficult 
to establish the level of the continuum beneath the PAH fea- 
tures, especially in the presence of 9.7 fim silicate absorption. 
The method presented here provides a systematic means of 
estimating this continuum — and therefore the strength of the 
PAH features — without the need for the observer to define the 
continuum level by hand. We model the flux density per unit 
frequency interval of the PAHs component as 



/■PAHs _ ft 

Jv -Jo 



PAHs 



r 



PAHs 



/PAHs(6.2/im) 



(34) 



where /^'^"^ = /^^'^(6.2y^m)-/^°"'(6.2/xm) is the estimated 
peak flux density of the 6.2 jim PAH feature obtained from 
the local continuum- subtracted observed spectrum, and 



f 



PAHs . 



E 



7?/i.o 



(A/A[,-Ai,/A) +7f 



(35) 



is the PAH emission template shown in Figure 8 constructed 
by summing over a series of Drude profiles with the param- 
eters presented in Table 1 . As suggested by the definition of 
the overall strength of the PAH template is held fixed 
during the continuum fit, but the strengths of the individual 
features are allowed to vary in the subsequent PAH fit (see 
§2.2 and the following paragraph). 

The central wavelengths and widths of our Ao > 3.3 [ivs\ fea- 
tures are adapted from those used in the analysis of Smith 
et al. (2007), and their relative strengths are derived using 







^0 


7r 


7^-' 


^z=0.25 


^.=0.5 


i 


Line 


(/iin) 


(10-2) 


(10-2) 


(10-2) 


(10-2) 


1 


H2S(8) 


5.053 




6.0 


6.0 


12.1 


2 


[Fell] 


5.340 


6.0 


6.0 


6.0 


12.1 


3 


H2S(7) 


5.511 


6.0 


6.0 


6.0 


12.1 


4 


[Fell] 


5.674 


6.0 


6.0 


6.0 


12. 


5 


H2S(6) 


6.109 


6.0 


6.0 


12.1 


12.1 


6 


[Fe 11] 


6.721 


6.0 


6.0 


12.1 


12.1 


7 


H2S(5) 


6.910 


6.0 


12.1 


12.1 


12.1 


8 


[Am] 


6.985 


6.0 


12.1 


12.1 


12.1 


9 


[Ne VI] 


7.652 


12.1 


12.1 


12.1 


12.1 


10 


H2S(4) 


8.025 


12.1 


12.1 


12.1 


12.1 


11 


[Arm] 


8.991 


12.1 


12.1 


12.1 


12.1 


12 


H2S(3) 


9.665 


12.1 


12.1 


12.1 


16.9 


13 


[SIV] 


10.511 


12.1 


12.1 


12.1 


16.9 


14 


H2S(2) 


12.279 


12.1 


12.1 


16.9 


16.9 


15 


[Ne II] 


12.814 


12.1 


12.1 


16.9 


16.9 


16 


[Nev] 


14.322 


16.9 


16.9 


16.9 


33.9 


17 


[Ne III] 


15.555 


16.9 


16.9 


22.3 


33.9 


18 


H2S(1) 


17.035 


16.9 


16.9 


33.9 


33.9 


19 


IS III] 


18.713 


16.9 


33.9 


33.9 


33.9 


20 


[Nev] 


24.317 


33.9 


33.9 


33.9 


33.9 


21 


[Oiv] 


25.890 


33.9 


33.9 


33.9 




22 


H,S(0) 


28.219 


33.9 


33.9 


33.9 




23 


[Slll] 


33.481 


33.9 


33.9 






24 


[Sin] 


34.815 


33.9 








25 


[Neiii] 


36.014 


33.9 









Note. — The FWHM = Aq7; of a particular atomic or molecular 
line is a function of the IRS module in which it is detected, and therefore 

depends on the redshift of the source. We give example values of 7,- for 
all lines at z = 0, 0.1, 0.25 and 0.5. Values of 7, are not given for lines 
that fall out of the IRS wavelength range at a given redshift. 

Pahfit'" to fit the low obscuration (rg 7 ss 0.24) mean star- 
burst galaxy spectrum created from 13 IRS spectra in Brandl 
et al. (2006). The 3.3 /im PAH feature central wavelength and 
width are obtained from Li & Draine (2001), and the peak flux 
density of this feature is derived from the ratio of the strengths 
of the PAH features at 3.3 and 6.2//m in the L-band and 
IRS spectra of the star-forming ULIRG IRAS 12112H-0305 
presented in Imanishi et al. (2006) and Armus et al. (2006). 
We include this 3.3 /zm PAH feature so that its contribu- 
tion to the integrated L-band (or IRAC 3.6 [im channel — see 
Fig. 8) flux density is explicitly included in the fitting method. 
As described in §2.2, after fitting the continuum, we fit the 
continuum-subtracted observed spectrum to refine our PAHs 
component (the properties of which are held fixed during the 
continuum fit). In this fit, the central wavelengths and FWHM 
of the PAH features are fixed to the values in Table 1, and the 
peaks are allowed to vary freely. This fit provides accurate 
luminosities for all PAH complexes, as weU as error estimates 
derived from the full covariance matrix. 

7. ATOMIC AND MOLECULAR LINE EMISSION 

In order to accurately fit the continuum and PAH features in 
an IRS spectrum, the contributions from atomic and molecu- 
lar lines must be properly fit as weU. Of course, the strength of 
these lines is scientifically interesting, but we caution against 
using our decomposition method to measure them since the 
fitted continuum around each line is not necessarily as ac- 
curate as one estimated by hand (a statement which cannot 
be made concerning the underlying continuum beneath PAH 
features). Therefore, although we obtain integrated flux val- 
ues for the lines we fit, we do not recommend using these 

See http://tiirtle.as.aiizona.edu/jdsmith/pahflt.php. 
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TABLE 3 
Dusty Galaxy Properties 



TABLE 4 
Supplementary Photometry 



Dl 



Galaxy 


z 


(Mpc) 


Class 


Reference 


NGC7714 


0.009 


39 


Starburst 


1,2 


NGC2623 


0.018 


78 


Obscured Starburst 


2 


PG 0804+761 


0.130 


610 


Quasar 


5 


Mrk463 


0.050 


222 


Seyfert 2 


4 


NGC6240 


0.024 


105 


Starburst+AGN 


3 


Mrk 1014 


0.163 


781 


Quasar+Starburst 


4 



REFERENCES. — [1] Brandl et al. (2004); [2] Brandl et al. (2006); [3] 
Armus et al. (2006); |4| Armus et al. (2004); |5| Hao et al. (2005a). 

Note. — The IRS spectra of the first three sources are from Spilzer pro- 
gram 14, and the spectra of the last three are from Spilzer program 105. 
Luminosity distances are calculated assuming a flat ACDM cosmology with 
Ho = 70 km s ' Mpc Qm = 0.3, and Qa = 0.7. 

measurements for scientific purposes (at least not without 
carefully checking the validity of the fitted local continuum 
around each Une). The one caveat to this recommendation 
where we do suggest using the fitted integrated line flux is for 
the [Ne II] 12.81 fim line, which our decomposition method 
is capable of deblending from the 12.7 /um PAH complex in 
low-resolution IRS spectra. 

We create our atomic and molecular emission line spectrum 
by estimating the peak flux density of each line in Table 2 
from the local continuum-subtracted IRS spectrum. The flux 
density per unit frequency interval of the lines component is 
modeled as a sum over unresolved Gaussian profiles having 
central wavelengths Aq, widths 7,, and peak flux densities f^fi 



/""" = EAoexp 



A- a;, 



Ai,7,■/(2^/2E2) 



(36) 



Since the Gaussian profiles are unresolved, the width of each 
line is a function of the IRS module in which it is detected (see 
Table 2 for values of 7 as a function of source redshift). As 
described in §6 for the PAHs component, after fitting the con- 
tinuum we perform a fit to the continuum-subtracted observed 
spectrum in order to refine the lines component (the parame- 
ters for which are held fixed in the continuum fit). In this fit, 
the central wavelengths and FWHM of the atomic and molec- 
ular emission fines are fixed to the values in Table 2, and their 
peak flux densities are allowed to range freely. From this fit, 
we obtain Une luminosities and their uncertainties — although, 
as noted above, we recommend measuring aU unblended line 
fluxes by hand. 

8. CONSTRUCTING DUSTY GALAXY SEDs 

8.1. Sample Selection 

Table 3 presents a summary of the galaxies in our sam- 
ple. They have been selected with the goal of demonstrat- 
ing the capabilities of our decomposition method on a variety 
of dusty galaxy SEDs, representative of the diverse proper- 
ties of the group. We include examples of unextinguished 
and obscured starbursts, a quasar, a Seyfert 2 galaxy, and two 
composite systems powered by both AGN and starburst ac- 
tivity (see §10.1 for details about each source). The first four 
galaxies in this sample serve as templates to aid in understand- 
ing the properties of the components from which composite 
systems are built. Our sample includes star-forming galaxies 
covering a range of infrared luminosities — from a minimum 
LiR « 5 X IO'^Lq for NGC7714 to more than a factor ~10 
higher for NGC 6240. While it is true that this sample does 



Galaxy 


Ultraviolet 


Optical 


Near-IR 


Far-IR 


Sub-mm 


NGC 7714 


1 


5 


5 


12,13 


13,17 


NGC 2623 


2 


6 


9 


12 


17,18 


PG 0804+761 


3 


3 


10,11 


14,15 




Mrk 463 


1 


7 


10 


14,16 




NGC 6240 


4 


8 


9 


12,16 


16,18 


Mrk 1014 


3 


7 


9 


14,15 





Rei ere.nCES. — [1| Kinney et al. (1993); [2] GALEX (2006); [3] Baskin & 
Laor (2005); [4] Sinith et al. (1992); [5] Lanfon et al. (2001); [6] Taylor et al. 
(2005); [7] Surace & Sanders (2000); [8] de Vaucouleurs et al. (1991); [9] Scoville 
et al. (2000); [10] 2MASS PSC; [11] Neugebauer et al. (1987); [12] Sanders et al. 
(2003); 113] Krugel et al. (1998); [14] Moshir & et al. (1990); [15] Haas et al. 
(2003); [16] Klaasetal. (2001); [17] Dunne et al. (2000); [18] Benford (1999). 

not include examples of 'normal' star-forming galaxies, our 
decomposition method is nonetheless vahd for these sources 
as well. In fact, our method is applicable to all galaxies — 
including AGNs, starbursts, normal, dwarfs, and ellipticals — 
so long as their SEDs are composed of emission from stars, 
dust, PAHs, atomic and molecular lines, and possibly an AGN 
accretion disk (i.e. our decomposition components). 

8.2. IRS Spectroscopy 

The IRS spectra of the sources in our sample were first pre- 
sented in the papers referred to in Table 3. All spectra in this 
paper have been extracted using the method described in Ar- 
mus et al. (2006). In brief, the IRS pipeline at the Spitzer Sci- 
ence Center was used to reduce the data, and one-dimensional 
spectra were extracted using the SMART package (Higdon 
et al. 2004). After scaling and stitching the individual orders 
together, each spectrum was scaled to match the MIPS 24 pm 
flux density if available (see §8.4), or the IRAS 25 /xm flux 
density otherwise. 

8.3. Supplementary Photometry 

The spectral coverage of the IRS low-resolution modules 
ends at an observed-frame wavelength of ~^38/im. Emission 
from typical cool component dust (f ~ 80 K) therefore only 
contributes to the last few microns, while emission from typ- 
ical cold component dust (T ^ 35 K) contributes negligibly 
at IRS wavelengths. To better constrain these cooler compo- 
nents which dominate the dust mass and frequently provide 
^50% of the total dust luminosity, we supplement our IRS 
spectra with f ar-IR to millimeter wavelength photometry from 
the hterature. Furthermore, the spectral coverage of the IRS 
modules begins at an observed-frame wavelength of ~5.2 /im. 
As seen in Figure 2, our source components all radiate sig- 
nificantly at shorter wavelengths. Thus, to better constrain 
these components, we supplement our SEDs with ultraviolet 
to near-IR photometry from the literature. Table 4 provides 
references to the literature for all photometry used in our de- 
compositions. 

In each wavelength range, we construct our SEDs using 
data which most closely samples the spatial region contained 
within the IRS slits. We treat photometric points that sam- 
ple larger spatial regions (i.e. for wavelength ranges where 
no data at a similar spatial resolution is available in the lit- 
erature) as upper-hmits to the flux density of the area cov- 
ered by the IRS spectrum. To include such a hmit as input 
to the fitting routine (which requires both the flux den- 
sity of each data point and its associated error), we assign the 
data to have an upper-Umit flux density /J: = and an uncer- 
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TABLE 5 

IRAC AND MIPS Nuclear Photometry 



IRAC MIPS 







3.6 


4.5 


24 


70 


160 


Galaxy 


PID 


(mJy) 


(mJy) 


(Jy) 


(Jy) 


(Jy) 


NGC7714 


59 


17.0 


13.5 


1.96 


6.67 


4.64 


NGC 2623 


32 


11.8 


13.5 


1.33 


16.4 


7.56 


PG 0804+761 


49 






0.19 


0.11 


0.033 


Mrk463 


32 


115 


156 


1.50 


1.90 


0.89 


NGC 6240 


32,3672 


33.7 


46.0 


2.71 


14.7 


9.08 


Mrk 1014 


32 


19.2 


23.2 








NOTE. —IRAC 


flux density 


errors are 


estimated to be 5% 


for Mrk 1 


014 and 



Mrk 463, and 10% for tiie other three sources. MIPS flux density errors are driven 
by the post-BCD calibration uncertainties and are ~ 10% for the 24 fim band and 
~20% for the 70 and 160 /nm bands. No color corrections are applied since they 
are much smaller than the uncertainties. The FID are the Spitzer program numbers 
for the IRAC and MIPS observations. 

tainty = /^^'^/3. If in addition to an upper-Umit to the flux 

density contained in the IRS sHt, we also have an estimate 
of the lower-limit (e.g. as estimated using a very small aper- 
ture measurement), we assign the data point to have a flux 

density of /,} = + /)])/2 with an associated uncertainty of 

fTj = (/J -/i)/3, so that the fit is constrained to fall between 
these two limiting values. 

8.4. IRAC and MIPS Photometry 

AU of the source components and the hot dust compo- 
nent radiate significantly from 3-5 /xm — i.e. between the JHK 
bands and the onset of the IRS spectral range (see Figs. 3 and 
7). To better constrain the fits in this critical region, we de- 
rive observed-frame 3.6 and 4.5 /im nuclear flux densities of 
the central point-source for five of six galaxies in our sample 
using archival IRAC images. This photometry is presented in 
Table 5 and is derived using the post-BCD products provided 
by the Spitzer Science Center. Flux densities for NGC 7714, 
NGC 2623, and NGC 6240 are calculated using a 5" radius 
circular aperture to match both the near-IR literature photom- 
etry and the IRS sUt-width. We do not apply any aperture 
corrections for these sources since they are all extended at 
these wavelengths. Mrk 1014 and Mrk 463 are both point-like 
at these wavelengths, so we calculate their flux densities using 
a 12.2" radius aperture, which requires no aperture correction 
to recover the full point-source value. Uncertainties are esti- 
mated to be 5% for Mrk 1014 and Mrk 463, and 10% for the 
other three sources. 

The IRAS and ISO beam-sizes are both much larger than the 
5-11" IRS slit-widths, so that far-IR photometry of nearby 
extended sources contains emission from regions outside the 
slits. We therefore derive mid-IR and far-IR nuclear flux den- 
sities of the central point-sources for five of six galaxies in our 
sample using archival MIPS images. This photometry is pre- 
sented in Table 5 and is derived using the post-BCD products 
provided by the Spitzer Science Center. We use Mrk 463 to 
derive aperture corrections since it is well-detected in all three 
bands and appears very point-like. To calculate the nuclear 
flux densities, we measure the emission from each source 
within small apertures (3", 8", and 20" radii at 24, 70, and 
160 /Ltm), and scale these values by the appropriate aperture 
correction (3.37, 4.09, and 2.95 for 24, 70, and 160 /xm). The 
S/Ji'teer beam-sizes at 24, 70, and 160 /im are ~6", 17", and 
39". The 24 /xm beam is therefore very similar to the 5-11" 
IRS sUts, justifying our choice to scale the IRS spectrum to 



this MIPS point. The beams at 70 /im and 160/im are signif- 
icantly smaller than the > 100" IRAS and ISO beams, so that 
the MIPS data provides a much better estimate of the nuclear 
flux density. Flux density errors are dominated by the post- 
BCD calibration uncertainties and are estimated to be 10% at 
24 /xm and 20% at 70 and 160^m. No color corrections are 
applied to the IRAC or MIPS photometry since aU corrections 
are much smaller than the uncertainties. 

9. DECOMPOSING DUSTY GALAXIES 

9.1. Spectral Decomposition Results 

In Figures 9 and 10 we present our decompositions of the 
nuclear SEDs of the dusty galaxies in our sample. In these 
fits, the MIPS photometric points are used to constrain the 
nuclear far-IR SEDs while the globally integrated IRAS, ISO, 
and sub-mm data provide upper-hmits to this nuclear emis- 
sion. The entire ultraviolet to millimeter wavelength SEDs 
are presented in Figure 9, and enlarged views focusing on the 
IRS regions are shown in Figure 10. In Figure 9, we also 
show the global fitted flux densities obtained using the full 
values of the large-beam IRAS, ISO, and sub-imn data. The 
temperatures of the cold and cool dust components obtained 
in the global flts are held flxed to these values in the nuclear 
decompositions. Note that the ISRF component contributes 
minimally (only in Mrk 463) since the small-aperture near-lR 
photometry used to match the IRS slits is dominated by emis- 
sion from young steUar populations and not evolved stars. 

9.2. Decomposition Parameters and Constraints 

The decomposition parameters obtained from the fits to the 
nuclear SEDs of the dusty galaxies in our sample are pre- 
sented in Tables 6 and 7. There are a total of 17 free parame- 
ters. The SB source component has two free parameters: dsB 
and fy^ determine the fractional SB contribution to the to- 
tal luminosity (see eq. [1]) and the 5500A optical depth (see 
eq. [6]). The ISRF source component has only a single free 
parameter: aisRF (see eq. [5]). The AGN source component 
has a total of two free parameters: Sagn and eagn, where 
the latter determines the fraction of the AGN component cov- 
ered by the obscuring screen (see eq. [9]). The optical depth 
through the obscuring clouds to the AGN component is fixed 
to f^*^*^ = 25. The hot and warm dust components each have 
three free parameters: a,-, fg^, and f/,- determine the compo- 
nent luminosities, 9.7 /xm optical depths, and dust tempera- 
tures via Ti = T{Ui) (see eq. [28]). The cool and cold compo- 
nents each have two free parameters: a,- and f/,- (see eq. [30]). 
Finally, there are two additional free parameters that control 
the water-ice and HAC contributions to the total opacity: ^ice 
and ?7hac (see eqs. [23] and [24]). 

We impose the following constraints to provide sensible 
limits over which the free parameters may range: (1) The 
magnitude of the radiation field energy density illuminating 
each dust component is constrained to give characteristic dust 
temperatures satisfying 20 K < Tioid < 50 K, 50 K < %ooi < 
150K, 150K < Twarm < 500 K, and 500K < TLot < 1500 K; 
(2) The maximum water-ice and HAC contributions to the to- 
tal opacity are constrained by their strengths in sources with 
very clean absorption spectra (see §4.8); (3) The extinction to 
the AGN source component is fixed to ti^ = 25 ; (4) The extinc- 
tion to the SB source component must satisfy ry < 5; (5) The 
optical depths through the screens obscuring the hot and warm 
components must not exceed the critical values beyond which 
their extinction-corrected luminosities would exceed the total 
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1 10 100 1000 1 10 100 1000 

X,^^^ (jlm) 1,^^, (jlm) 

Fig. 9. — Spectral decompositions of tlie ultraviolet to millimeter wavelength SEDs (black solid lines and filled circles) of the nuclear regions of the dusty 
galaxies NGC 7714, NGC 2623, PG 0804+761, Mrk 463, NGC 6240, and Mrk 1014. Arrows indicate upper or lower-limits (see §8.3). Unfilled circles indicate the 
position of the MIPS 24 /im point to which the IRS spectra are scaled. IRAC and MIPS photometric points ai'e circled. The various decomposition components are 
shown as indicated in the figure legend. Additionally, the total fitted SEDs obtained using the globally integrated far-IR data (dashed 'Total' line) are displayed. 
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Fig. 10. — Closeup view of the IRS wavelength ranges for the nuclear decompositions presented in Figure 9. The various decomposition components are 
shown as indicated in the figure legend. Also shown ai'e the l-cr uncertainties in the IRS spectra {dark shaded regions) and the formal l-cr error contours for the 
total fit (light shaded regions). Emission lines and the wavelength ranges over which silicate emission and absorption is observed are also labeled. 
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TABLE 6 

Nuclear SED Decomposition Parameters and Characteristic Dust Temperatures 



Galaxy 


/doP 


dof 


T b 
-/cold 


T J 


7"- b 
J- warm 


T b 


^warmc 
'9.7 


'9.7 


'V 


~ e 

^agn 


'/ice 


f 

'/HAL. 


NGC7714 


1.33 


382 


31 


71 


165 




0.19 




0.97 
















±1 


±1 


±1 




±0.03 




±0.02 








NGC2623 


2.46 


383 


29 


53 


206 




3.52 




4.81 




0.1 


0.1 








±1 


±1 


±1 




±0.06 




±0.11 








PG 0804+761 


2.77 


363 


42 


150 


400 


1410 


0.78 


1.22 





0.25 


0.050 


0.044 








±2 


±6 


±3 


±10 


±0.01 


±0.03 




±0.01 


±0.012 


±0.035 


Mrk463 


1.90 


364 


30 


55 


211 


1130 


1.30 


2.06 


1.12 


1 


0.012 


0.062 








±2 


±2 


±1 


±6 


±0.05 


±0.01 


±0.26 




±0.004 


±0.009 


NGC 6240 


1.77 


369 


29 


61 


193 


1260 


3.64 





4.68 


1 


0.067 


0.1 








±1 


±1 


±1 


±30 


±0.07 




±0.08 




±0.007 




Mrk 1014 


2.00 


362 


31 


71 


195 


1500 


0.14 


2.70 


5 


0.49 


0.031 


0.1 








±1 


±1 


±1 




±0.02 


±0.04 




±0.01 


±0.005 





Note. — Formal statistical parameter uncertainties are given beneath eacli parameter value. If no uncertainty is given, the corresponding parameter 
was pegged at a limiting value. The 5500 A optical depth to the AGN component is fixed to f,')'^^ - 25 for all fits. 

^ Total weighted reduced value (see eqs. [2] and [3]) for the indicated number of degrees of freedom (dof).'^ Characteristic dust temperature (see 
§4.3) determined from the fitted value of the magnitude of the illuminating radiation field energy density, t/, .'^ 9.7 fim optical depth through the screen 
obscuring the warm and hoi dust components (see eqs. [28]).*^ F-band (5500 A) optical depth through the screen obscuring the SB source component 
(see eqs. [6] and [7]).'^ Fraction of the AGN source component covered by the obscuring screen (see eq. [9]).^ Ratio of the 6.1 jum water-ice and 
6.85 ^m HAC optical depths to the 9.7 ^m optical depth of each dust component (see eqs. [23] and [24]). 



dust luminosity; (6) The ratio of the unextinguished luminos- 
ity of the SB component to the luminosity of the PAHs compo- 
nent must be between the values obtained for our unobscured 
and obscured template starbursts NGC 7714 and NGC 2623; 
and (7) The hot dust covering fraction (i.e. the ratio of the un- 
extinguished hot dust luminosity to the unextinguished AGN 
luminosity) must be > 0.5, consistent with observations and 
corresponding to torus opening angles > 30°. 

9.3. Derived Quantities 
9.3.1. Source and Dust Component Luminosities 

Apparent and extinction-corrected source and dust compo- 
nent luminosities from the nuclear decompositions are pre- 
sented in Table 7. The apparent luminosity of each component 

is calculated from 



Li = 4ttD 



/^source 



Li G ^source 5 



(37) 



where Di is the luminosity distance to the galaxy (see Ta- 
ble 3), and ^source and L^ust are the sums of the luminosities 
of all source and dust components, respectively. Extinction- 
corrected luminosities are calculated from 



L\ = Ai,Dl 



f 



Li ^ ijs, 



(ly = ) '^i ^source '-'i ^ ^.^source) 
T,(i^) 1 Q^j'i'dust l-i € I-dust. 



(38) 



where the T, are the extinction terms for each component (i.e. 

Tagn = [ 1 - eagn] + cagn exp[-TAGN] for the AGN component 
and T, = exp[-T,] for the warm, hot, and SB components). 

9.3.2. Dust Component Masses 

The mass of dust in each of the cold, cool, warm, and hot 
nuclear decomposition components are provided in Table 8. 
We generalize the expression used to calculate dust masses for 
a single grain-size and composition (e.g. see eq. [4] of Klaas 
et al. 2001) to the case in which grains are distributed both in 
size and composition. The mass of dust in each component is 
therefore given by 



M,=Z) 



/i(Ao) 



£UAo)e-^.<^«) 



E / (39) 
^ J 3 riH da 



which is the product of the number of hydrogen nucleons 
required to power the fitted dust component and the total 
dust mass per hydrogen nucleon for the adopted dust model 
(summed over the graphite and silicate compositions). In this 
equation, Ao is an arbitrary wavelength at which the terms are 
evaluated, and the densities of graphite and silicate are taken 
to be pgra = 2.24 and psii = 3.50 g"^ cm~^. For comparison to 
the dust masses and luminosities obtained from the nuclear 
decompositions, we present the corresponding properties for 
the global decompositions in Table 9. 

9.3.3. Comparison With Previous Results 

Klaas et al. (2001) fit the far-IR to sub-nmi wavelength 
SEDs of a sample of ULIRGs using dust-modified blackbod- 
ies with temperatures typically ranging between 28 and 50 K. 
Included among their sources are two galaxies which are also 
in our sample, Mrk 463 and NGC 6240, for which they fit sin- 
gle (3 = 2 dust components with temperatures of 40 and 33 K, 
respectively. In a similar study of quasars presented in Haas 
et al. (2003), the far-lR SED of PG 0804+761 is fitted with a 
47 K dust component. The temperatures we derive from our 
decompositions are all consistent with these values. Our de- 
composition of Mrk 463 includes contributions from 30 K and 
55 K dust components, while our fit to NGC 6240 includes 
emission from dust at 29 K and 61 K. Finally, our decomposi- 
tion of PG 0804+761 is dominated in the far-lR by emission 
from 42 K dust. For all three sources, the temperatures of the 
single component fits in the literature fall between the tem- 
peratures of the dust in our multi-component fits. Klaas et al. 
(2001) also fit the A > 60 /im SED of NGC 6240 with multi- 
ple (3 = 2 dust components and obtain temperatures of 3 1 and 
41 K (they also include a 10 K component to fit to sub-mm 
upper- hmits). Again, these temperatures are quite comparable 
to ours, noting that our warmer component is at a somewhat 
higher temperature in order to fit to the shorter wavelength 
data included in our fits. 

As illustrated by these decompositions, multi-component 
fits including dust at different temperatures are rarely unique 
since several components can usually be added together to 
approximate the SED of emission from a smaller number of 
components. In fact, dust in galaxies is likely to be distributed 
throughout many regions with various heating conditions, and 
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TABLE 7 

Nuclear SED Decomposition Source and Dust Component Luminosities 



Source 



Dust 



Galaxy 


Source 

(IQIOLq) 


"ISRF 


"SB 


"SB 


Q^AGN 


"^AGN 


^dust 

(lO'^L©) 


"cold 


cicool 


dwarm 


5' 


Shot 


"hot 


NGC7714 


0.78 





1 


3.8 








3.25 


0.344 


0.394 


0.262 


0.275 










±0.03 




±0.05 


±0.1 






±0.04 


±0.011 


±0.006 


±0.004 


±0.004 






NGC 2623 


0.72 





1 


27 








22.5 


0.136 


0.792 


0.072 


0.194 










±0.05 




±0.09 


±2 






±0.2 


±0.010 


±0.010 


±0.001 


±0.004 






PG 0804+761 


99 











1 


1.27 


62.5 


0.036 


0.146 


0.457 


0.639 


0.361 


1 




±2 








±0.03 


±0.04 


±1.1 


±0.002 


±0.003 


±0.011 


±0.013 


±0.014 


±0.02 


Mrk463 


8.2 


0.12 


0.23 


1.03 


0.65 


18.1 


74 


0.068 


0.165 


0.470 


0.703 


0.297 


1 




±0.6 


±0.03 


±0.05 


±0.08 


±0.05 


±1.4 


±2 


±0.003 


±0.005 


±0.027 


±0.038 


±0.009 


±0.03 


NGC 6240 


2.61 





>0.99 


24.8 


<0.01 


0.29 


45.2 


0.257 


0.528 


0.198 


0.514 


0.017 


0.017 




±0.07 




±0.04 


±0.7 




±0.01 


±0.7 


±0.011 


±0.009 


±0.007 


±0.020 


±0.001 


±0.001 


Mrk 1014 


120 





0.018 


0.49 


0.98 


1.86 


337 




0.583 


0.336 


0.352 


0.081 


0.66 




±3 




±0.001 


±0.01 


±0.03 


±0.04 


±2 




±0.005 


±0.003 


±0.003 


±0.002 


±0.02 


Note. — Apparent and extinction-corrected luminosities of each 


source and dust component 


, in the nuclear decompositions are given as fractions of Lsource 


and Ldust (see eqs. [37] 



and [38]). Formal statistical uncertainties for the fitted luminosities are given beneath their values. The nuclear dust luminosity of Mrk 1014 is uncertain by ' 
between its nuclear and global values) and its cold component luminosity is undetermined since no MIPS data is available to constrain its far-IR SED. 



'20% (i.e. the difference 



TABLE 8 

Nuclear SED Decomposition Dust Masses 









Mwai-m 


Mho, 


Galaxy 


{UfM.J 


(lO'Af .J 


(lO-Af,, 1 


{M . ) 


NGC 7714 


4.26 ±0.13 


4.90 ±0.04 


5.55 ±0.05 





NGC 2623 


19.8 ±1.4 


335 ±2 


9.3 ±0.2 





PG0804+761 


1.53 ±0.07 


0.91 ±0.01 


2.95 ±0.04 


1.75 ±0.09 


Mrk 463 


23.7 ±0.9 


185 ±3 


101 ±13 


7.9 ±0.2 


NGC 6240 


67.3 ±2.6 


201 ±1 


69±6 


0.049 ±0.003 


Mrk 1014 




736 ±5 


333 ±3 


1.13 ±0.03 



Note. — Formal statistical uncertainties for the fitted masses of the dust components 
in the nuclear decompositions are given along with their values. The mass of cold dust 
in the nucleus of Mrk 1014 is tmdetermined since no MIPS data is available to constrain 
its far-IR SED. 



is therefore distributed over a continuum of temperatures. 
Any multi-component decomposition therefore serves only to 
characterize these distributions of temperatures. Additionally, 
if either or both of the restrictions that dust is optically-thin 
to its own emission and that its far-IR emissivity index is 
characterized by /3 = 2 imposed above are lifted, yet differ- 
ent dust temperatures may be obtained for the same SEDs. 
For example, Klaas et al. (2001) perform fits in which the 
dust optical depth and emissivity are free parameters, and ob- 
tain much higher temperatures — 52 and 57 K for Mrk 463 and 
NGC 6240, respectively. Given the multiple factors determin- 
ing the temperatures derived from spectral decompositions 
(e.g. fixed versus variable (3), we emphasize that caution must 
be taken when comparing temperatures derived from different 
methods to ensure that they are in fact comparable. 

The total far-IR luminosity from dust inferred from a multi- 
component decomposition is largely independent of the num- 
ber of components and the choice of their characteristic tem- 
peratures (since the luminosity is just the integral under the 
SED). The total mass of dust, however, depends upon these 
choices (since the mass is a strong function of temperature). 
Since the total mass of dust in galaxies is dominated by low- 
temperature material, decompositions with multiple compo- 
nents (including very cold dust) produce much higher dust 
masses. For example, Klaas et al. (2001) estimate the mass 
of dust in Mrk 463 to be greater than 5 x 10* Mq and less 
than 1.2 X 10^ Mq as derived from their single and multi- 
component fits, respectively. Similarly, they estimate that 



the mass of dust in NGC 6240 is between 3.6 x 10^ and 
5.8 X IO^Mq. Our decompositions yield total globally inte- 
grated dust masses of 2.8 x lO^M© and 1.6 x lO^M© for the 
two sources, relatively consistent with the previously derived 
masses (vvilhin the uncertainties of the methods). Klaas et al. 
(2001) estimate the H2 mass of NGC 6240 to be 3.7 x 10'°Mq 
(assuming M[H2] = 4.6Lco). This value, along with our esti- 
mate of the total dust mass in NGC 6240, gives a gas-to-dust 
ratio of ~230, in reasonable agreement with the canonical 
value of ~165 from Li (2005). 

9.3.4. Diagnostic Ratios 

In Table 10, we present several 'diagnostic' ratios for each 
galaxy derived from the nuclear decomposition masses and 
luminosities presented in Tables 7 and 8. The first entry in 
this table gives the ratio of the total dust luminosity to the to- 
tal gas mass. The total gas mass is calculated from the fitted 
value of the total dust mass and the gas-to-dust ratio of our 
adopted dust model: Mgas/Mdust = l ^Me/Mdust = 124 (Li & 
Draine 2001). The definition of this quantity is motivated by 
the fact that two sources with equivalent values of Ldust or Lir 
may have completely different physical properties (e.g. a very 
massive but inactive system may have the same dust luminos- 
ity as a relatively small but highly active galaxy). The total 
dust luminosity is therefore a degenerate quantity which does 
not alone provide great insight into the nature of a galaxy. 
Normalizing the dust luminosity by the dust mass, however, 
does provide a means of breaking this degeneracy. This ratio 
characterizes the radiative efficiency of a galaxy, with smaller 
Ldust/ Mgas ratios indicating lower-efficiencies and cooler dust 
temperatures and larger ratios indicating higher-efficiencies 
and warmer dust temperatures. 

Also included in Table 10 are the ratios of the extinction- 
corrected SB and AGN source component luminosities to the 
total nuclear dust luminosity. A ratio of unity or greater in 
these quantities indicates that the total observed dust emis- 
sion could be heated by the fitted source component. A ratio 
significantly greater than unity indicates that either the source 
extinction correction was over-predicted, or that the dust cov- 
ering fraction is less than unity. In the latter case, the fraction 
of the source which is covered by dust is obtained from the 
inverse of the source-to-dust ratio. A ratio significantly less 
than unity indicates that either the extinction to the source 
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TABLE 9 

Global SED Decomposition Dust Properties 





^dust 






^cold 


^cool 


Galaxy 


(IQiO Lq) 


"cold 


Q^cool 


(10« Mq) 


(lO^Mo) 


NGC7714 


4.50 ±0.05 


0.530 ±0.010 


0.280 ±0.006 


9.1±0.1 


4.8±0.1 


NGC 2623 


30.3 ±0.9 


0.365 ± 0.036 


0.581 ±0.019 


72± 13 


331±3 


PG 0804+761 


62.6±3.8 


0.050 ± 0.004 


0.140 ±0.0 10 


2.2 ±0.1 


0.88 ±0.03 


Mrk463 


74.5 ± 1.3 


0.076 ±0.014 


0.153 ±0.007 


27 ±5 


173 ±8 


NGC 6240 


61.2±3.7 


0.461 ±0.039 


0.381 ±0.040 


164 ±21 


197 ±60 


Mrkl014 


414 ±23 


0.188±0.015 


0.469 ± 0.054 


303 ± 18 


730 ±280 



Note. — Global SED decomposition properties derived from fitting the total far-IR emission (see 
Fig. 9). All quantities are analogous to those in Tables 7 and 8. The temperatures of the cold and cool 
components are the same as given in Table 6. Formal statistical uncertainties for the fitted properties are 
given along with their values. 



TABLE 10 

Nuclear SED Decomposition Diagnostic Ratios 





Dust 


Starburst 


AGN 


PAHs 


Galaxy 


^dust/^gas^ 




^AGN/^'gls'^ 
^AGN/^dust'^ (LqMq) 


ipAHs/A^gas'' 
ipAHs/idust'' (^-0*^0) 



NGC 7714 


60.4 ±2.0 


0.912 ±0.018 


55.3 ±1.9 






0.052 ±0.002 


3.17±0.13 


NGC 2623 


78.2 ±4.7 


0.877 ±0.028 


68.5 ±4.6 






0.021 ±0.001 


1.65±0.11 


PG 0804+761 


3270 ±160 






2.00 ±0.06 


6530 ±327 






Mrk463 


233 ±10 


0.1 14 ±0.003 


26.5 ±0.9 


2.00 ±0.06 


466 ±16 






NGC 6240 


52.4 ±2.1 


1.43 ±0.02 


75.1 ±2.8 


0.017 ±0.001 


0.88 ±0.03 


0.032 ±0.002 


1.65±0.11 


Mrk 1014 




0.173 ±0.001 




0.661 ±0.001 




0.012 ±0.001 





Note. — Formal statistical uncertainties for the diagnostic ratios are given along with their values. Mass-normalized values are not given for Mrk 1014 since no 
MIPS data is unavailable to constrain its far-IR SED (and therefore Mgas). This absence of data also results in a ~20% uncertainty for the value of Ldust for Mrk 1014 
(i.e. the difference between its nuclear and global values), and therefore in all quantities normalized by this total dust luminosity. 

^ Ratio of the total integrated dust luminosity to the total gas mass.^ Ratios of the extinction-corrected SB source component luminosity to the total dust luminosity 
and total gas mass.*^ Ratios of the extinction-corrected AGN source component luminosity to the total dust luminosity and total gas mass."* Ratios of the total PAHs 



component luminosity (integrated over all features) to the total dust luminosity and total gas mass. 



component is underpredicted, or that additional unmodeled 
source emission is present which has not been accounted for 
in the decomposition (e.g. the apparent near-lR luminosity of 
a deeply embedded source may be small in comparison to a 
relatively unextinguished source, even though the two may be 
intrinsically similar in strength). In addition to the source-to- 
dust luminosity ratios, we also calculate source-to-dust mass 
ratios, where the extinction-corrected SB and AGN luminosi- 
ties are given as a fraction of the total gas mass. This quantity 
indicates the relative 'dustiness' (in terms of mass) of a source 
compared to others of the same luminosity. 

9.3.5. PAH Feature Luminosities 

The total apparent PAH luminosity of each source (inte- 
grated over all features), normalized to both the total dust lu- 
minosity and the total gas mass, are presented in Table 10. 
Additionally, the fitted ratios of the luminosities of the four 
primary PAH complexes (see Table 1) to the total PAH lumi- 
nosity are presented in Table 1 1 . For reference, we also pro- 
vide the corresponding ratios for the mean starburst galaxy 
spectrum from Brandl et al. (2006) used to create our PAH 
template, and the median ratios for the star-forming SINGS 
galaxies from Smith et al. (2007). With the exception of 
Mrk 1014, the fitted PAH luminosity ratios of the four sources 
in our sample with well-detected PAH emission are fairly con- 
sistent with the values of the mean starburst and SINGS tem- 
plates (see §10.1 for more on the PAH emission of Mrkl014). 

The apparent ratios of PAH emission from a galaxy may 
differ from those of the low-obscuration templates in Table 1 1 
for several reasons. First, the size of a PAH (i.e. the number 



of carbon atoms) and its ionization state can alter the rela- 
tive strengths of its emission features (see Draine & Li 2001). 
Additionally, the apparent ratios may vary as a result of ex- 
tinction from intervening dust. In Table 12, we present the 
ratios of the 6.2, 7.7, and 11.3 /xm PAH complexes to the total 
PAH luminosity for the mean starburst galaxy from Table 1 1 
after passing through screens of dust with various values of 
Ti)j. In addition to the feature ratios, we also provide correc- 
tion factors used to convert the total apparent PAH luminosity 
into the actual emitted PAH luminosity for a particular value 
of Tgj estimated from the feature ratios. 

10. DISCUSSION 

10.1. Detailed Descriptions of Decompositions 

NGC 77 14 — Our low-obscuration template starburst is a 
massive barred spiral galaxy that makes up the western com- 
ponent of the interacting system Arp284. The compact 
ultraviolet-luminous starburst nucleus of NGC 77 14 has very 
low extinction, and there is no evidence at any wavelength 
for the presence of an AGN. The IRS spectrum of this galaxy 
(see Brandl et al. 2004) exhibits a rising continuum, strong 
PAH feature and line emission, and no perceptible 9.7 /zm sil- 
icate absorption. As suggested by the relatively fiat contin- 
uum near 10 /im, all parameters from our fitting are consistent 
with this starburst galaxy being relatively unobscured at all 
wavelengths. The warm dust and SB source components have 
fitted values of rg j w 0.19 and ry = 0.97, respectively. Puxley 
& Brand (1994) estimate a screen extinction to the ionized gas 
within NGC 7714ofAvR:!0.86±0.13 using optical and near- 
IR recombination fines, a value which is entirely consistent 



Decomposing Dusty Galaxies. I. 



19 



TABLE 1 1 
PAH Feature Luminosities 



Galaxy 


^6.2/ipAHs 


^T-v/^-PAHs 


-f-ll.s/^PAHs 


^Iv/^-PAHs 


NGC7714 


0.114 ±0.002 


0.396 ± 0.004 


0.1 14 ±0.002 


0.067 ± 0.004 


NGC 2623 


0.088 ±0.002 


0.407 ±0.005 


0.074 ±0.001 


0.026 ±0.002 


PG 0804+761 










Mrk463 










NGC 6240 


0.1 15 ±0.005 


0.371 ±0.009 


0.105 ±0.002 


0.069 ± 0.004 


Mrk 1014 


0.112±0.013 


0.211 ±0.014 


0.067 ±0.008 


0.051 ±0.016 


Mean Starburst 


0.10 


0.39 


0.11 


0.08 


Median SINGS 


0.11 


0.42 


0.12 


0.06 



Note. — Ratios of the luminosities of each of the four main PAH complexes to the total PAH 
luminosity. Formal statistical uncertainties for the fitted properties are given along with their values. 
Also shown are the ratios for the mean starburst galaxy spectrum used to derive our PAH template 
(see Table I) and the median values for the SINGS sample of star-forming galaxies presented in 
Smith et al. (2007). 



TABLE 12 

Extinguished Mean Starburst PAH Feature Luminosities 



79.7 


^6.2/i'PAHs^ 




i-ll.s/i-PAHs^ 


^PAHs' ^PAHs 


0.25 


0.101 


0.402 


0.106 


1.09 


0.5 


0.104 


0.411 


0.099 


1.19 


1 


0.111 


0.429 


0.087 


1.40 


2 


0.123 


0.462 


0.067 


1.91 



* Complex-to-total PAH luminosity ratio for the indicated screen 
T9.7.'' Unextinguished-to-apparent PAH luminosity ratio for the indicated 
screen r^j. 



with our SB component optical depth. For our adopted dust 
model, Ay « 0.86 corresponds to rgj « 0.07-0.22, depending 
on the influence of scattering in the total effective opacity (the 
latter value is for zero scattering contribution, as expected for 
spherical symmetry). This range is again completely consis- 
tent with the extinction we find for the warm dust component. 

The extinction-corrected SB-to-dust ratio given in Table 10 
is Ljg / Ldust = 0.91, indicating that the fitted source component 
is capable of powering nearly all (i.e. 91%) of the total dust 
emission. Although we attempt to calculate the true far-IR 
dust luminosity of the nuclear point-source, we may nonethe- 
less overestimate this quantity since our MIPS photometry 
may still contain flux from outside the nucleus. Such an over- 
estimate could give rise to the fitted slight source-luminosity 
deficit (i.e. there is not enough source emission to explain the 
observed dust emission). Given this possible explanation, and 
the fact that the deficit is small (< 10%), we conclude that the 
observed source-to-dust ratio of NGC 7714 is consistent with 
the total system being powered by star-formation. 

Our decomposition of this galaxy therefore indicates that 
the warmest dust in pure low-obscuration starbursts typically 
has a temperature ~165K, and that emission from ~1000K 
dust is therefore negligible. Also, the observed PAH luminos- 
ity in NGC 7714 is ~5.2% of its total nuclear dust luminosity, 
providing an estimate of the typical total PAH strength in low- 
obscuration starbursts. Smith et al. (2007) find a median PAH- 
to-dust luminosity ratio of ~10% for the 'normal' (i.e. non- 
starburst) star-forming galaxies in the SINGS sample, with 
the ratio ranging between 3.2 and 16% for the sample. Thus, 
the ratio measured for NGC 77 14 is consistent with this range, 
with the nearly factor of two difference from the median value 
likely indicating real differences in the emission from galaxies 
with different levels of star-formation. Note also that the PAH 
ratio L11.3/L7.7 0.28 for NGC 7714 is very similar to the 
value of this ratio for the unobscured mean starburst template 



(~0.27). The PAH emission in NGC 7714 therefore appears 
to be relatively unobscured, consistent with the results from 
other extinction indicators. 

The total global dust emission in our sample (i.e. obtained 
using IRAS and ISO far-IR data) ranges from ~1.0 to 1.4 
times the nuclear emission (i.e. obtained using MIPS data). 
Langon et al. (2001) estimate that only 10-50% of the stellar 
emission in NGC 7714 emerges from the 2" x 2" (^0.37pc^) 
nucleus, and therefore infer that 50-90% of the far-IR IRAS 
emission likely emerges from the extended disk of the galaxy. 
They further reason that the strength of the extra-nuclear far- 
IR emission is likely to be towards the lower end of this range, 
since the ultraviolet observations used to quantify the stellar 
emission strength are prone to missing emission from extin- 
guished populations. In comparing the results of our global 
and nuclear spectral decompositions of NGC 7714, we calcu- 
late that at least ~50% of the cold dust component luminosity 
(and 30% of the total dust emission) emerges from extended 
regions, consistent with the Langon et al. (2001) prediction. 

NGC 2623 — Our extinguished template starburst galaxy is 
a late stage merger with long tidal tails and a r'/'* nuclear 
stellar Ught profile (Wright et al. 1990; Scoville et al. 2000). 
The radio continuum (Condon et al. 1991) and mid-IR (Soifer 
et al. 2001) emission is similarly dominated by a compact nu- 
clear source. The IRS spectrum of NGC 2623 is dominated by 
deep silicate absorption features at 9.7 and 18 /im and a very 
steep continuum beyond 20 fj.m. The fitted optical depths to 
the warm dust and SB source components of NGC 2623 are 
T9,7 w 3.5 and Ty » 4.8. Note that Brandl et al. (2006) find 
an apparent 9.7 fim optical depth of ~1.5, which is more than 
a factor two lower than our fitted value. This difference oc- 
curs as a result of the method they employ to determine the 
apparent optical depth — ^namely interpolating a smooth con- 
tinuum above the silicate absorption feature and using the ra- 
tio of this 'unextinguished' continuum to the observed value 
to infer T9 7. In contrast, the unextinguished continuum in the 
dust model used in our decomposition method is not smooth 
at 9.7 (Um, but instead features a silicate emission feature (the 
presence of which is consistent with observations — e.g. see 
the fit for PG0804H-761), therefore requiring a higher optical 
depth to obtain the same extinguished continuum. 

We find no evidence for the presence of an AGN contri- 
bution to the SED of NGC 2623, consistent with the simi- 
lar conclusion in Risaliti et al. (2000) based upon hard X- 
ray observations. The extinction-corrected 5B-to-dust ratio 
given in Table 10 is Lsg/Ldust = 0.88, indicating that the fitted 
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source component is capable of powering most of the total 
dust emission. As with NGC 7714, it is likely that this ratio is 
indicative of an entirely star-formation powered system (tak- 
ing into account the likely small overprediction of the nuclear 
dust luminosity due to the MIPS resolution). Like NGC 7714, 
the mid-IR spectrum of NGC 2623 also exhibits significant 
emission from PAHs, albeit a factor ^^2.5 weaker. The PAH 
ratio L11.3/L7.7 ~ 0.18 for NGC 2623 is also lower than the 
value for NGC 77 14, with the lower ratio being consistent 
with PAH emission obscured by a screen of dust with rg j ~ 1 
(as determined by comparison with the obscured mean star- 
burst template properties in Table 12). For this optical depth, 
the apparent PAH luminosity of the mean starburst template 
is only ~71% its true emitted value. Using this correction 
factor, we estimate that the true PAH-to-dust luminosity ra- 
tio of NGC 2623 is ipAHsZ-Ldust ~ 0.029, a factor ~1.8 lower 
than NGC 7714. See §10.2 for a detailed comparison of the 
properties of the low obscuration system NGC 7714 and the 
obscured starburst NGC 2623. 

PG0804+761— Our Type-1 template AGN is a variable 
source displaying ^40% near-lR and 10.6 /im brightness fluc- 
tuations over the course of a decade (Neugebauer & Matthews 
1999). The first deep infrared observations of this source were 
obtained with ISO (Haas et al. 2003), which revealed a power- 
law spectral shape and a very warm infrared SED peaking at 
^30 //m. Analysis of its IRS spectrum by Hao et al. (2005a) 
revealed the presence of strong 9.7 and 18 /im amorphous sil- 
icate emission features. In our decomposition, these emission 
features are produced by grains in the warm dust component. 
If we assume that the warm component is directly illuminated 
by the PG 0804+761 accretion disk (Lagn = 1.26 x IO'^Lq), 
then equation (12) with Ucui, L^/Airr^ indicates that this 
dust is located ~19pc from the nucleus. Jaffe et al. (2004) 
analyze 10 /xm VLTl interferometric observations of the nu- 
cleus of the nearby Type-2 AGN NGC 1068, and conclude 
that its torus emission must be confined to a region < 2pc in 
size. Thus, some of the optically-thin warm component emis- 
sion from PG 0804+761 may emerge from regions beyond the 
torus, perhaps from clouds in the narrow-hne-region. Note 
that this distance to the warm dust is actually an upper-limit to 
the true distance, since the emission may also originate from 
indirectly illuminated torus clouds (see Fig. 6) at r < 19pc. 

The unextinguished AGA^-to-dust ratio in the PG 0804+761 
decomposition is ~2 (see Table 10), indicating that the AGN 
source component is more than capable of heating the total 
observed dust emission. The inverse of this ratio provides an 
estimate of the total dust covering factor (i.e. the fraction of 
sky covered by the torus as seen from the nucleus), which is 
/torus ~ 0.5 (note that this value is at the constraint imposed 
during the fit — see §9.2). This covering factor is consistent 
with a torus opening angle of ^30° (i.e. sinf?torus = /torus, 
where 6'torus is the angle extending from the mid-plane to the 
top of the torus as viewed from the accretion disk). Based 
upon the statistics of Seyfert galaxies (i.e. the relative num- 
bers of Type-1 and 2 galaxies), Schmitt et al. (2001) find ap- 
proximately twice as many Type-2 AGN, and hence estimate 
6'torus = 45°. More recently, Hao et al. (2005b) study the AGNs 
in the Sloan Digital Sky Survey and conclude that (?torus ~ 30°, 
consistent with the value obtained from our decomposition. 

Jaffe et al. (2004) find that the 8-13jum VLTI spectrum of 
NGC 1068 can be decomposed into a T > 800 K dust compo- 
nent obscured by a T9.7 = 2. 1 ± 0.5 screen and af = 320±30K 
component behind a rg j = 0.3 ± 0.2 screen. The temperatures 



and screen optical depths of these components are comparable 
to those obtained in our spectral decomposition (see also the 
properties of Mrk463). PG 0804+761 has the worst reduced 
of the galaxies in our sample, due primarily to the poor 
fit to its IRS spectrum between 7.5-10 fjm (although the fit to 
the rest of the SED is quite good). Jaffe et al. (2004) report 
a similar problem in this wavelength range in their fit to the 
spectrum of NGC 1068, and conclude that the fit is improved 
using high-temperature calcium aluminium silicate dust in- 
stead of the standard olivine-type siUcates (the opacity of the 
9.7 /im feature of the former species begins near 9/im as op- 
posed to the ^8 /im onset for olivine). We finally note that our 
decomposition of the ultraviolet and optical emission is con- 
sistent with a largely uncovered AGN accretion disk — with 
only ~25% covered — consistent with a Type-1 source. 

Mrk463 — Our Type-2 template AGN is a merging system 
with two nuclei separated by ^4" (Mazzarella et al. 1991). 
The portion of the IRS spectrum obtained with the short- 
low module is likely dominated by Mrk463e (since the slit 
of this module is only ~3.6" wide and the observation tar- 
geted the eastern component), while the portion obtained with 
the wider-sht of the long-low module likely contains emission 
from both nuclei (although the optically-bright western com- 
ponent may not contribute significantly at these wavelengths). 
Mrk463 has a luminous steep-spectrum radio core, and broad 
lines are seen in scattered optical (Miller & Goodrich 1990) 
and direct near-lR light (Goodrich et al. 1994; Veilleux et al. 
1997). Like the quasar PG 0804+761, the mid-lR spectrum 
of this galaxy shows no discernible PAH emission features 
(Armus et al. 2004). Silicate absorption at 9.7 /xm is clearly 
seen in the spectrum of Mrk463, producing a mid-lR contin- 
uum dominated by a broad emission-like feature at ~8 /xm, 
characteristic of self-absorbed siUcate dust. In our decompo- 
sition, this feature is primarily produced by the warm com- 
ponent which has a screen optical depth rg.v « 1.3. Like 
the quasar PG 0804+761, Mrk463 has a moderately obscured 
(r9 7 w 2.1) hot dust component. The hot component opti- 
cal depth and the deviations of the fit to the IRS spectrum 
from 7.5-9 /xm are both similar to the reported properties of 
PG 0804+761 and NGC 1068 described above. UnUke the fit 
to the quasar PG 0804+761, the AGN component of Mrk463 
is completely covered by obscuring clouds — i.e. eACN = 1 — as 
is expected for a Type-2 source. 

The AGA^-to-dust ratio for Mrk463 is 1.93 (similar to the 
value for PG 0804+761), consistent with the AGA'^ source 
component powering the total observed dust emission, and 
implying a torus opening angle 6'torus ~ 3 1 ° . There is evidence 
for ISRF and SB contributions in the decomposition, although 
these are bolometrically weak compared to the power of the 
AGN « 5%). Although Mrk463 and PG 0804+761 have sim- 
ilar AGA^-to-dust ratios, they have very different quantities of 
cooler dust. Mrk463 has 16 times more cold dust by mass 
(see Table 8), and it also has a factor 14 times smaller dust 
luminosity-to-mass ratio (also indicative of cooler mean dust 
temperatures — see Table 10). Similarly, the ratio of the unex- 
tinguished AGA^ component luminosity to the total gas mass 
is a factor 14 smaller for Mrk463 (i.e. it has 14 times more 
gas mass for the same accretion disk luminosity). In the stan- 
dard AGN unification scenario, the only difference between 
a Type-1 source (such as PG 0804+761) and a Type-2 source 
(such as Mrk463) is the orientation of the obscuring torus 
with respect to the observer. Thus, within this picture, Type-1 
and Type-2 AGNs should have similar cold dust masses rela- 
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tive to their total masses, and they should therefore have sim- 
ilar AGN luminosity-to-mass ratios (since the accretion disk 
luminosity and the total mass of dust do not depend on ori- 
entation). In §10.4, we suggest that the relative orientations 
of the accretion and host galaxy disks may explain why this 
simple consequence of the unification scenario does not hold 
for some galaxies. We therefore emphasize that the SEDs of 
individual AGNs may appear falsely inconsistent with the uni- 
fication scenario due to local geometrical factors, even if uni- 
fication holds for AGNs as a class. 

NGC6240 — The first of our two composite sources is an 
interacting system containing two nuclei separated by ~1" 
(so that both nuclei are contained within all IRS slits). With 
an LiR w 6.5 x lO^L© NGC6240 is technically a LIRG, al- 
though it is often treated as a ULIRG since it shares most 
properties with other members of the class. Based upon its 
optical nuclear spectrum, NGC 6240 is classified as a LINER 
(Armus et al. 1989), and X-ray observations (Komossa et al. 
2003) provide clear evidence for the presence of a pair of 
AGNs located behind significant columns of absorbing neu- 
tral gas (A^H = 1-2 X 10^"*cm"^). The mid-lR spectrum of 
NGC 6240 is extremely rich, displaying strong PAH features, 
both low and high ionization lines (e.g. [Ne II] 12.8/xm and 
[Ne V] 14.3 ^m) as well as strong emission from molecu- 
lar hydrogen (see Armus et al. 2006). Like NGC 2623, the 
mid-lR continuum of this galaxy is shaped by strong absorp- 
tion from silicate grains at 9.7 and IS/xm, as indicated by 
the warm component screen optical depth of T9.7 « 3.6. The 
near-IR continuum of NGC 6240 is well-fitted by a combi- 
nation of apparently unobscured weak hot dust and extin- 
guished (tv = 4.7) SB component emission. The detection 
of [Ne V] 14.3 /im in the high-resolution IRS spectrum of 
NGC 6240 presented in Armus et al. (2006) and the evidence 
for X-rays are both consistent with the presence of hot dust in 
our decomposition — indicating that a (perhaps small) fraction 
of the near-lR emission is powered by an AGN. 

Based upon Chandra X-ray observations. Lira et al. (2002) 
conclude that the AGN in NGC 6240 likely contributes be- 
tween 30-50% of the bolometric luminosity. However, as re- 
ported in Armus et al. (2006), the small [Ne V]/[Ne n] and 
[Ne V] /LiR flux ratios are both consistent with an AGN con- 
tribution of only 3-5% of the bolometric luminosity. Sim- 
ilarly, the apparent hot dust emission in our decomposition 
makes up just ~2% of the total dust luminosity (compared to 
> 36 and 30% for PG0804H-761 and Mrk463, respectively), 
so that any contribution to the observed near- and mid-lR from 
an AGN is very small. On the other hand, the 5B-to-dust ra- 
tio for NGC 6240 is 1.43, indicating that the total observed 
dust emission could be powered by the observed SB compo- 
nent. This galaxy therefore presents quite a puzzle, whereby 
its appearance changes radically depending upon where you 
look — although X-rays suggest powerful AGNs, data at op- 
tical and longer wavelengths do not require such a presence. 
We note that even though the hot component in the decompo- 
sition is unobscured, it could still be associated with a deeply 
obscured AGN in a clumpy geometry where most of the hot 
dust is too obscured to be seen, and we view only the most 
unobscured portions (see §10.3). If this scenario is true, the 
actual hot dust contribution to the total bolometric luminosity 
could easily be pushed higher. 

MrklOM — Our second composite source is a radio-quiet 
infrared-luminous QSO with broad optical emission lines 



(FWHM H/3 > 4000 km s"'). Mrk 1014 displays twin tidal 
tails, indicative of a recent interaction and merger (MacKenty 
& Stockton 1984). It is a relatively warm far-IR source, with 
the peak of its SED located around 70 /um. Its mid-IR spec- 
trum is characterized by a nearly power-law continuum, with 
no obvious silicate emission or absorption features, and weak 
PAH emission (see Armus et al. 2004). This power-law con- 
tinuum is rather remarkable (assuming the continuum is ther- 
mal in origin) since a nearly featureless SED must be con- 
structed from emission components which are known to have 
significant features (i.e. silicates). The decomposition which 
provides a solution to this puzzle is consistent with Mrk 1014 
being a composite source, containing a combination of dust 
components seen in both our template starbursts and AGNs. 
Like the AGNs, Mrk 1014 contains a moderately obscured 
(T9.7 « 2.7) hot dust component and a much less obscured 
warm component (tq 7 « 0. 14). Like the starbursts (and unlike 
the quasar), Mrk 1014 also contains bolometrically significant 
emission from cooler dust, comprising ^60% of its total dust 
luminosity (we argue in §10.4 that some of this far-IR emis- 
sion is likely powered by the AGN). 

The extinction-corrected AGA^-to-dust and 5B-to-dust lumi- 
nosity ratios of Mrk 1014 are 0.66 and 0.17, respectively, in- 
dicating that the AGN is capable of powering ~66% of the 
observed dust emission from the galaxy. As fitted, the AGN 
and SB components are together capable of powering 83% of 
the observed dust emission, thereby requiring the presence of 
at least an additional 17% source luminosity to power the re- 
maining dust emission (note also that the total nuclear dust lu- 
minosity from Mrk 1014 is uncertain by ~20% since no MIPS 
data is available to constrain its far-IR SED). If all of this un- 
detected source luminosity emerges from the AGN accretion 
disk (e.g. if we underpredict the AGN component luminosity 
as a result of too rigid a constraint on the /lof-to-AGA^ luminos- 
ity ratio), the AGN would account for 83% of the bolometric 
luminosity of Mrk 1014. If, on the other hand, the remaining 
undetected source luminosity is entirely produced by a star- 
burst (e.g. if we underpredict the SB luminosity as a result 
of missing emission from highly embedded regions), then the 
AGN would account for 66% of the bolometric luminosity. 

Armus et al. (2007) create diagnostic diagrams based upon 
mid-IR spectral lines (see their Figs. 5-8) which are all con- 
sistent with an AGN fraction between 50-90% for Mrk 1014. 
Additionally, Boiler et al. (2002) conclude that Mrk 1014 is 
dominated by an AGN and not star-formation based upon X- 
ray observations. If we take the PAffi'-to-dust luminosity frac- 
tions of the galaxies NGC 2623 and NGC 6240 to be represen- 
tative of the range of values in obscured starbursts, then the 
PAHs-to-d\xst ratio of Mrk 1014 suggests that between ^38 
and 57% of its bolometric luminosity is powered by obscured 
star-formation (assuming all PAH emission comes from star- 
formation). Given our conclusion above that the starburst con- 
tribution is between 17 and 34%, it is therefore likely either 
that we are underestimating the total starburst contribution 
based upon the source components (which seems improbable 
given the corroborating X-ray and emission line evidence), or 
that some of the observed PAH emission may actually be pow- 
ered by the AGN (see §10.4). We note as well that the PAH 
feature luminosity ratios for Mrk 1014 (see Table 1 1) — in par- 
ticular the 7.7 iim complex — differ significantly from the tem- 
plate starbursts, perhaps indicative of an AGN contribution 
giving rise to non-standard PAH emission. 

10.2. Comparison of Starburst SEDs 
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From our modest sample of two starburst galaxy decom- 
positions, it is apparent that starburst SEDs do not form a 
homogeneous group. Comparing the fits of NGC 2623 and 
NGC 7714, noticeable differences include a 30% higher max- 
imum dust temperature, a factor 20 times higher maximum 
optical depth, a factor 2.5 times weaker PAHs-to-dust lu- 
minosity ratio, and comparatively weaker [Ne III] 15.56 /im, 
[S III] 18.71 A«m, [S III] 33.48 ^m, and [Si II] 34.82/xm line 
emission. If the nuclear emission from each starburst emerges 
from distinct star-forming clouds consisting of embedded 
stars surrounded by cocoons of dust (see Fig. 1), then many 
of the observed variations in the SEDs of NGC 7714 and 
NGC 2623 can be understood as resulting from: (1) differ- 
ences in the total number of star-forming clouds (i.e. if indi- 
vidual cloud complexes in each galaxy have similar intrinsic 
luminosities, then the total luminosity of each galaxy is de- 
termined by the total number of clouds it contains), and (2) 
the mean optical depth through an individual cloud. There 
are, of course, many other factors shaping the SEDs of star- 
bursts as well — e.g. age and metalhcity. The simple model 
described in this section is therefore not intended to provide 
a comprehensive explanation of all starburst properties, but 
rather a structure within which to understand several general 
properties. 

If we assume for simphcity that each cloud in a starburst 
has constant gas density, then the optical depth through a 
cloud is T(r) = Tcioud('-/'-cioud), where Tdoud and rdoud are the 
total optical depth and radius of the cloud. As described in 
§5.1, direct emission from stars heats a shell of dust out to a 
radius r^heu, corresponding to the optical depth Tgj", where 
the cloud becomes optically thick to the heating stellar pho- 
tons (i.e. essentially the outer edge of the photodissociation 
region). Using the equation for the optical depth above, we 

nhtaln r2623 /,J714 _ /_7714 / 2623 w 2623 /^7714x tt 

oDtam rj,^eii/' shell - V cloud/ ^cioudA'^cioud/'^cioudi- "^rc, wc as- 
sume that T9 j " is the same for all clouds since it depends only 
on the properties of dust and the illuminating radiation field. 
With Tcioud Twann (i-C- using the warm component optical 
depth as a proxy for the total cloud optical depth), this gives 
ddi /d™ ~ (0.05)(4^„?/r™ ). If the clouds in each galaxy 
are similar in size, this suggests that ''^heii / ''shin < ^ (i-^- '^^e 
warm dust shell from Fig. 6 is geometrically more extended in 
NGC 7714 than NGC 2623), so that stellar photons penetrate 
to much larger radii in the star-forming clouds of NGC 7714. 
Furthermore, the fractional cloud volume containing dust at 
the temperature of the warm component should be smaller 
in NGC 2623 than NGC 77 14, and consequently the mass of 
warm dust relative to the total mass should be smaller as 
well. Indeed, from Table 8 we find that (M„arm/Afdust)77i4 ~ 
1.3 X 10"^ and (Mwann/^^dust)2623 ~ 4 X 10"^, Consistent with 
this prediction. 

For r < Tsheii, energy conservation dictates that the temper- 
atures of grains within a cloud scale roughly as T oc r"'/-^. 
This, coupled with the fact that r^^ < r^^^^, impUes that 
the mean temperature of such grains should be higher in 
NGC 2623 than NGC 7714— a prediction which is consistent 
with its 25% higher warm component temperature. The lu- 
minosities per unit mass of dust at the temperatures of the 
warm components in the two galaxies scale approximately 
as (r^arm/T'wrm)'* ~ (1.25)'' « 2.4 (i.e. a mass of warm dust 
from NGC 2623 will be ^2.4 times more luminous than an 
equivalent mass from NGC 7714). This ratio, along with the 
absolute warm component dust masses for the two galaxies, 
suggests that the ratio of the warm component luminosities 



should be L^^^^ /L^^"* W (f2623 /fn\Mrj^2623 /y^^7714^ 4 
aiiuuiu uc i^^y^rm / ^warm V ^ warm / ^ warm / v"warm/ ^"warmV ~ ^■ 

The actual fitted ratio is 4.7, in reasonable agreement with this 
prediction. Furthermore, the Ldust/^H ratio from Table 10 
is ~30% higher in NGC 2623 than NGC 77 14. If the warm 
dust in NGC 7714 were at the temperature of the warm dust 
in NGC 2623, the total dust luminosity of NGC 7714 would 
increase by ^40%, and the resulting dust luminosity-to-mass 
ratios of the two sources would differ by only ~7%. Thus, 
the different values of Ldust/^gas are directly related to the 
higher warm dust temperature in NGC 2623, which itself is 
Ukely directly related to the higher optical depth through the 
star-forming clouds in that galaxy. 

The apparent PA//i-to-dust luminosity ratio is a factor ^2.5 
times lower in NGC 2623 than NGC 7714. We estimate that 
extinction to the PAHs component in NGC 2623 reduces its 
luminosity by a factor 1.4 (see §10.1), so that the unextin- 
guished PAHs-to-dust luminosity ratios for the two galaxies 
differ by a factor 1.8. With this correction, we also find that 
i-PAHs/A^gas is only a factor 1.4 greater in NGC 7714, as op- 
posed to the observed factor 1.9. As described above, the 
warm-to-total dust mass ratio is less in NGC 2623, so that the 
volume within which PAH molecules are heated is smaller 
for this source. Therefore, some portion of the reduced PAHs 
luminosity giving rise to the differences in ratios described 
above likely results from such a difference in heating geome- 
tries. Using the ratio of the extinction-corrected LpAHs/Afgas 
values derived above, we estimate that this geometrical effect 
must decrease the observed PAHs luminosity of NGC 2623 by 
a factor 1 .4 (i.e. assuming that the full difference is caused by 
this geometrical effect). Note that this factor is similar to the 

value [(M„arm/Mdust)2623/(M„aim/Mdust)7714](i™??^/iwal4) ~ 

1.4, where the first term accounts for the difference in vol- 
ume within which PAHs are heated, and the second term ac- 
counts for the difference in heating intensity. With these two 
correction factors of 1.4, the PAHs-to-dust luminosity ratio 
of NGC 2623 differs by only ~30% from that of NGC 7714. 
Recalling that the the dust luminosity of NGC 7714 would in- 
crease by ~40% if its warm dust were at the same temperature 
of the dust in NGC 2623, the fi4i/s-to-dust luminosity ratios 
of the two galaxies differ by less than 10%. 

The differences in the apparent PAH emission from our two 
starburst templates may therefore be understood as resulting 
from a combination of: (1) extinction to the PAHs, (2) the 
geometry of the PAH emitting regions, and (3) the method 
of normalizing the PAHs component luminosity. We suggest 
that the ratio LpAHs/Mgas is a somewhat better tracer of the 
true strength of PAH emission than LpAHs/idust, since it is 
less affected by temperature effects. This latter ratio, or a 
related quantity such as L6,2/^fir, is often used to quantify 
the amount of star-formation in galaxies (e.g. Peeters et al. 
2004). Given that both our template starbursts are believed 
to be entirely powered by star-formation, it is clear that there 
is significant variation in the relative strength of PAH emis- 
sion in pure starbursts. We therefore caution that great care 
must be taken when using any of these metrics to derive ab- 
solute quantities of star-formation, since both extinction and 
geometrical effects may result in different ratios for pure star- 
burst galaxies. 

10.3. Hot Dust Emission from Clumpy AGN Tori 

If the obscuring dust surrounding an AGN accretion disk 
is arranged in a smooth (i.e. non-clumpy) toroidal distribu- 
tion, then the ~3-8/xm SEDs of Type-1 sources should be 
dominated by emission from relatively unobscured hot dust 
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Fig. 11. — Edge-on (upper-left) and face-on (upper-right) views of a 
clumpy torus comprised of a toroidal distribution of discrete clouds. Di- 
rectly illuminated cloud faces (i.e. the sources of hot emission) are obscured 
from most viewpoints by dust on the shaded side of each cloud. SEDs of the 
fitted NGC6240 and PG 0804+761 hot components are shown beneath the 
schematic torus model representing their geometry (solid lines), along with 
model SEDs derived assuming that a fraction of the hot component emission 
(50% for NGC6240 and 75% for PG 0804+761) is obscured by a tv = 25 
screen, with the remainder emerging unobscured (dashed lines). 



(since this hot dust is viewed directly by the observer, with 
no intervening material from the torus along the line-of-sight 
to obscure it). If, however, the obscuring structure is instead 
composed of discrete clouds distributed in a roughly toroidal 
shape (i.e. a clumpy torus), then hot dust emission is not ex- 
pected to emerge unobscured. As illustrated in Figure 1 1, the 
directly illuminated (and therefore hotter) faces of clouds in a 
clumpy torus always point radially towards the nucleus, and 
therefore away from the observer. Thus, emission from these 
hot cloud faces, and therefore hot component emission from 
such sources, is obscured by cooler dust on the shaded sides 
of clouds. Our decomposition of the quasar PG 0804-H761 in- 
cludes a moderately obscured (T9.7 « 1 .2) hot dust component, 
consistent with such a clumpy obscuring torus. 

The SED of the hot dust component from the PG 0804H-76 1 
decomposition is shown in the lower-right panel of Figure 1 1 . 
Also shown is a curve representing emission from dust at the 
temperature of the PG0804H-761 hot component, with 75% 
obscured by a = 25 screen (i.e. equivalent to the screen ob- 
scuring the AGN component), and the remaining 25% emerg- 
ing unobscured. While this model with non-uniform coverage 
does not exactly reproduce the SED from the fully covered 
model (in particular at shorter wavelengths), it nonetheless 
produces a very similar spectrum which would also provide 
an adequate fit to the PG0804H-761 SED (with suitable small 
adjustments to the AGN and warm components). Thus, the 
SED of hot component emission from a highly obscured ge- 
ometry with a few clear lines-of-sight may appear similar to 
the emission from a fully covered but less obscured geometry. 
For a traditional torus made from a smooth dust distribution, 
hot dust emission from a face-on source would emerge unob- 
scured with strong silicate emission features (since there are 
no shaded sides of clouds to obscure the hot emission). Our 
decomposition of PG0804-H761 does not show such strong 
features from the hot component, suggesting that its hot dust, 
and possibly that of most Type-1 sources, may therefore be 
significantly obscured by the shadowed sides of clumpy torus 
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Fig. 12. — If the disk of the host galaxy and the AGN accretion disk are 
co-planar (top), very little dust in the disk of the galaxy is illuminated by the 
AGN (i.e. the radiation cones from the accretion disk do not intersect the 
host galaxy disk). In contrast, if the host galaxy and AGN accretion disks are 
orthogonal (bottom), the AGN heats dust over a large area of the galaxy disk 
(i.e. the radiation cones do intersect the host galaxy disk), resulting in far-IR 
and possibly PAH emission. Broad arrows label the line-of-sights giving rise 
to the indicated properties of the emergent SEDs for the various geometries. 



clouds. 

The lower-left panel of Figure 11 depicts the fitted 
NGC 6240 hot dust component SED (unobscured in the de- 
composition), along with a model in which 50% of its hot dust 
is covered by a screen with Ty = 25 and the remainder emerges 
unobscured. As with PG 0804-1-761, this partially covered 
model does not precisely reproduce the fitted hot dust com- 
ponent. It is, however, similar enough to the fitted component 
to provide a reasonably good fit (i.e. since the hot component 
of NGC 6240 contributes in a region with strong contributions 
from several other components, the small deviations between 
the two models would not affect the decomposition greatly). 
The near equivalence of these two models demonstrates that 
emission from a highly obscured, but only partially covered, 
geometry may appear quite similar to emission from a fully 
covered geometry at much lower levels of obscuration. There- 
fore, despite the small apparent optical depth derived from the 
fit to NGC 6240, it is possible that the hot component in this 
source is actually much more significantly obscured (i.e. we 
may only see direct emission from a few clear lines-of-sight, 
while the majority of the hot dust is completely obscured). 

10.4. Origin of the Far-IR Emission in AGNs 

Since emission from dust at far-IR (A > 50/xm) wave- 
lengths is likely to emerge unobscured from galaxies, the av- 
erage properties of the far-IR SEDs of AGNs (e.g. their rela- 
tive cold dust luminosities and masses) should be similar for 
all sources (assuming their nuclei are similarly structured), in- 
dependent of their orientation-based (i.e. Type- 1 versus Type- 
2) classification. As described in §10.1, the two template 
AGNs in our sample have very different far-IR properties (e.g. 
Mrk463 has a dust luminosity-to-mass ratio ^15 times that of 
PG0804H-761), in sharp contrast to these expectations. Far- 
rah et al. (2003) use a two-component (starburstH-AGN) tem- 
plate decomposition to conclude that some optically classified 
AGNs with strong far-IR emission likely have significant con- 
tributions from starbursts. While this is likely true for many 
AGNs with cool far-IR SEDs, it does not appear to hold for 
Mrk463 (see §10.1). As illustrated in Figure 12, we instead 
suggest that the far-IR emission from some AGNs (including 
Mrk463) is powered by the AGN itself and originates from 
cool dust in the disk of the host galaxy. We note here that 
this suggestion is based upon the SEDs of only two AGNs, so 
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further study using a much larger sample (to be presented in 
a subsequent paper) is needed before conclusions are drawn 
about the origin of far-IR emission from the class of AGNs. 

While it is often assumed that an AGN accretion disk and 
the disk of the galaxy hosting the AGN are co-planar (as in the 
top of Fig. 12), this geometrical assumption is not necessarily 
always true. After the merger of two galaxies, the resulting 
accretion and host galaxy disks may eventually come to rest 
in the same plane (e.g. due to torques), but there is a period of 
time after the merger during which their relative orientation is 
essentially random. If the two disks are tilted with respect to 
each other, a significant quantity of dust in the disk of the host 
galaxy may be heated by the AGN out to large radial distances 
(see the bottom of Fig. 12). Assuming that the galaxy disk is 
not a significant source of opacity, a galaxy therefore appears 
as a Type- 1 or Type-2 AGN depending upon the relative ori- 
entation of the obscuring torus with respect to the observer 
(i.e. the standard AGN unification scenario). Additionally, 
depending on the relative orientation of the host galaxy and 
AGN accretion disks (but independent of the orientation of 
the observer), the source may exhibit weak or strong far-IR 
emission (i.e. "Warm" or "Cold" far-IR SEDs as labeled in 
Fig. 12). In this model, Mrk463 and PG 0804+761 are there- 
fore "Cold" and "Warm" far-IR sources, respectively, with 
emission arising from geometries similar to those in the bot- 
tom and top of Figure 12, respectively. We note that the ge- 
ometry suggested here has a similar impact on the SEDs of 
galaxies as the 'warped-disk' geometry from Sanders et al. 
(1989). 

As described in §10.1, we estimate that ~35% of the bolo- 
metric luminosity of Mrk 1014 is powered by star-formation. 
If all of this luminosity emerges in the form of cool dust, 
this still leaves ^30% of its 70K emission component to be 
powered by its AGN. We therefore suggest that the geom- 
etry of Mrk 1014 is hkely to be similar to that of Mrk 463 
(i.e. the bottom of Figure 12). Additionally, the luminosity 
of the Mrk 1014 PAHs component compared to its inferred 
starburst fraction suggests that between ^20-50% of its PAH 
emission is powered by its AGN (i.e. derived by compari- 
son with the properties of the obscured starbursts NGC 2623 
and NGC 6240— see §10.1). Schweitzer et al. (2006) argue 
against AGN powered far-IR emission in a sample of QSOs 
observed with Spitzer based largely on the fact that PAHs are 
destroyed in a hard AGN radiation field. However, in the ge- 
ometry of Figure 12, the radiation field gradually softens as 
the photons propagate through the galaxy disk, so that at some 
radius PAH molecules illuminated by the AGN may survive. 
Clearly, more work is required (e.g. radiative transfer calcu- 
lations with realistic PAH destruction models) to verify if this 
scenario is indeed reahstic, and we do not dispute that some 
far-IR emission in AGNs is powered by starbursts, but the 
suggested mechanism to generation far-IR and PAH emission 
in AGNs as a result of the interaction of the accretion disk 
with its host galaxy should be considered as well. Further- 
more, we warn against the blanket assumption that all PAH 
and far-IR emission is directly associated with star-formation, 
as this could result in an overestimation of the true level of 
star-formation activity in galaxies. 

11. SUMMARY 



We have presented a new multi-component spectral energy 
distribution decomposition method. Our principal goal for 
this paper is to demonstrate the effectiveness of the method 
for studying composite dusty galaxies, although we note that 
the method is quite general and may be appUcable to much 
broader selections of galaxies. To demonstrate the efficacy of 
the method, we have applied it to the ultraviolet through mil- 
limeter wavelength SEDs of the nuclei of the unobscured and 
obscured template starbursts NGC 7714 and NGC 2623, as 
well as the template Type-1 and Type-2 AGNs PG 0804+761 
and Mrk 463. We find that the pan-spectral SEDs of these 
sources are all consistent with being entirely powered by star- 
formation and accretion, respectively. We compare the SEDs 
of the two starburst galaxies and demonstrate that the dif- 
ferences in their properties (including the luminosities and 
masses of their dust components and the strength of their PAH 
emission) can largely be explained in terms of the mean opti- 
cal depth through the dusty cocoons surtounding their newly 
formed stars. In comparing the template AGNs, we find that 
the far-IR SEDs of the two sources differ significantly, which 
we suggest may result from variations in the relative orienta- 
tion of their host galaxy and AGN accretion disks. 

We also decompose the nuclear SEDs of the composite 
sources NGC 6240 and Mrk 1014. Our decomposition finds 
very little evidence for an AGN contribution in NGC 6240 
(~2%), with the total dust emission in that galaxy being pow- 
ered by a starburst with mean optical depth Ty w 5. We 
note, however, that we cannot rule out a larger contribu- 
tion from a deeply embedded AGN visible only in X-rays. 
The SED of Mrk 1014 is consistent with a galaxy which is 
~65% AGN powered, with a ~35% contribution from star- 
formation. Like Mrk 463, we find that Mrk 1014 has a far-IR 
excess (compared to the Type-1 AGN PG 0804+761), and at- 
tribute this to cold dust in the disk of its host galaxy which is 
heated directly by its AGN. We also estimate that up to 50% 
of the detected PAH emission in Mrk 1014 may be powered 
by radiation from the AGN that has softened as it propagates 
through the galaxy disk (and therefore does not destroy the 
fragile PAH molecules). In a future paper, we will apply this 
decomposition method to a much larger sample of starbursts, 
AGNs, LIRGs, and ULIRGs observed as part of the Spitzer 
GTO program in order to provide a suite of templates to aid 
in understanding both the local dusty galaxy population and 
the high-redshift LIRGs and ULIRGs now being uncovered 
by Spitzer. 
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APPENDIX 

WEIGHTING OF NON-UNIFORMLY SAMPLED DATA 

The sampling term, A, in equation (3), is required to produce sensible and fair fits of non-uniformly sampled data, as is the case 
when coarsely sampled photometry is joined together with spectral data. If such a term is excluded (i.e. = 1 for all A^), then 
will be dominated by the more finely sampled spectral data, simply because it contains many more data points. To establish a 
weighting scheme for non-uniformly sampled data, we define the quantity 

niRS = Mrs/ log(A^^fjA'JfiS„), (Al) 
which is the number of spectral samples per logarithmic wavelength interval in an IRS spectrum. Here, Nirs is the total number 
of spectral samples in the IRS spectrum, and Ajjjj^ = min(AiRs) and AjJ^^^ = max(AiRs), where Airs is the array of wavelength 
points in the IRS spectrum. We then divide the full SED into Nuns wavelength bins (typically 10), distributed in even logarithmic 
wavelength intervals, with each bin Z e [0, 1 , . . . ,Mbins - 1] containing minimum and maximum wavelengths 

/A \ /A \ 

Ajni„ = min(A/) = A^in I I and A^^^ = max(A;) = Amin I I , (A2) 

\ '^min / \ ^min / 

where Amin = min(A) and Amax = max(A) are the minimum and maximum wavelength in the SED having wavelength array A. 
Each bin is then assigned a total weight based upon the sampling of the IRS spectrum 

A' „ i„„ f ^'mnx \ Mrs l0g(Amax/ Amin) 

^^- = "^"^°H^j = ^log(A-/A-)- ^^3) 

The weight of bins containing no data points is redistributed to the nearest bins containing data above and below, with the weight 
going to each bin being a linear function of the proximity to the bin with no data; i.e. for a bin /' containing no data, the weight 
distributed to the bins /+ and l'_ above and below is 

^Ki^J" = Kins'^^ and 6Air'J-=A(i^JjrrjF- (^4) 

The weighting in each bin is then distributed to the data points within that bin, so that for a data point at wavelength k which is in 
bin / ^ 

Aie, = %^, (A5) 
TV; 

where Ni is the number of data points in bin I. Thus, the A^^ represent the number of effective spectral samples that each data 
point represents so that the sampling over the entire SED is equivalent to the IRS sampling. Finally, we normaUze and scale the 
weighting function 

A, = Mot^^, (A6) 

so that = Mot = Mpiiot+MiRs, the total number of data points in the spectral energy distribution, as would be the case if the 

weighting function were taken to be unity at all wavelengths (i.e. A^ = 1 for all Xk). 

IRS STATISTICAL FLUX DENSITY UNCERTAINTIES 

We derive the statistical flux density uncertainty at each wavelength element of an IRS spectrum by differencing two spectra 
obtained at different nod positions of each slit. The initial estimate of the uncertainty at each wavelength is 

'^*= iV2Y ' ^^'^ 

where f\ and /a are the observed flux densities in nod positions 1 and 2, and the two factors of \/2 both correct for the introduction 
of noise resulting from the subtraction, and account for the reduction in noise of the final flux densities resulting from averaging 
the two nods. At each wavelength, this array contains a single sample of the true uncertainty drawn from the probability distri- 
bution function having standard deviation (j{Xk). We estimate the true value of the statistical uncertainty at each wavelength by 
calculating the standard deviation of the points in a window of width W (typically around 20 wavelength bins, with the window 
size decreasing near the edges of the spectrum), and assigning this value as the standard deviation of the central wavelength, i.e. 

k+W/l 
i=k-Wl2 

where {5f) « is the mean value of 5/ for the points within the window. Since we are approximating a function which in 
principle should be sampled on an infinitely fine grid, our estimated a array will always show some residuals from our finite grid 
size. As a result, the a array does not always have the expected property that ^68% of the 5 fk are enveloped by the l-cr contour. 
We therefore perform one final step, smoothing our estimated a array with a window of width W, and scaling the entire array 
until it has the expected behavior that the l-a error contour envelops 68% of the 5fk. 
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